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(involving the recalculation of the whole of the arithmetical work) was pointed out by Mr Arthur Berry, 
who examined the paper in proof. I am deeply grateful to him for the discovery. Owing to Mr Young’s 


call to military service he has been unable to reconstruct the paper. 


Accordingly for the paper in its 


present form Miss Elderton and I are responsible, not only for the calculations, but for the much modified 
and considerably augmented analysis. K. P. 


The Torsion resulting from Flexure in Prisms with 
Cross-sections of Um-aaial Symmetry. 


INTRODUCTION 


The problem of the flexure of a beam under the action of a single load stands 
today in much the same state as Saint-Venant left it half a century ago, nearly 
all of the more recent work having been devoted to the problem of flexure with 
continuous systems of loading. There has not even been any notable additions to 
the forms of the cross-section of the beam discussed by Saint-Venant and the 
rectangular and elliptical cross-sections may be taken as the only cross-sections of 
simple mathematical form for which full solutions exist. It appears, for instance, that, 
although the general solution of the problem on Saint-Venant’s assumptions includes 
the case of a beam having a uni-symmetric cross-section, no such case has been 
discussed analytically, far less numerically. Yet the case of general asymmetry in the 
cross-section possesses peculiar interest and has recently become of real importance 
in connection with the blades of aeroplane propellers. With the high speed of 
rotation to which the aeroplane propeller attains there is a strong flexing force in 
the air resistance and, as we shall see, the flexure of a beam of asymmetrical 
cross-section produces a torsion of a by no means negligible amount. To the 
designer of the aeroplane propeller this is of the utmost importance for even a slight 
change in the ‘angle of attack’ of a propeller-blade might seriously alter the effi- 
ciency and speed of the machine. It may well be that more potent causes exist 
for the torsion of a propeller-blade subjected to a continuous system of pressure 
and not merely to a single end-load, but to gain some knowledge of the nature 
of the torsion produced by flexure in the simpler case of the present paper must 
be regarded as the first step towards the investigation of the more difficult problem. 
There is the further limitation to the generality of the discussion of the following 
pages that we shall treat only of cross-sections which have an axis of symmetry 
perpendicular to the plane of flexure. The problem of complete asymmetry still 
awaits investigation. 

The torsion spoken of above which arises in the flexure of a beam of asym- 
metrical cross-section can in general terms be seen to arise from two contributory 
causes. Firstly there is the distortion of the cross-sections of the beam which 
Saint-Venant showed to be'a necessary accompaniment of flexure. With an 


asymmetrical cross-section the distortion, it is clear, will be likewise asymmetrical. 
1—2 
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If then we fix the end-section in almost* any way which will allow the distortion 
of the theory to take place, this distortion and particularly that part of it which 
accompanies the curvature, called by Thomson and Tait anticlastic, will produce a 
rotation of the beam as a whole, together with what we may call an ‘impure’ 
torsion. The latter we shall term the torsion due to anticlastic curvature. It 
may be noted in passing that even in the case of a beam of symmetrical cross- 
section, there will arise a rotation of this nature if the end fixing be asymmetrical - 
in character. 


This torting effect of the distortion of the sections is however only one portion 
of the total torsion; with it must be combined the effect due to the distribution 
of shears{ in every cross-section of the beam—a distribution which is necessarily 
asymmetrical if the cross-section be asymmetrical. It is the main object of this paper 
to investigate this torsion due to shear with more particular reference to the case 
where the cross-section is a circular sector, complete or curtate. The solution, of 
the corresponding torsioh problem exists for this cross-section, and with numerical 
reductions such as we should expect in work done by Saint-Venant{. Arithmetical 
work is even more necessary in the flexure problem and no small amount of care 
has been devoted in this paper to make clear the analysis by numerical calculation, 
whenever that was possible. In addition, more perhaps with a view to strengthening 
the authors’ faith in the theoretical results than to proving them to the complete 
satisfaction of others, a series of rough experimental tests were carried out with 
various beams of circular sector cross-section. In the course of the discussion notes 
will be found of these observations. They are given for what they are worth and 
make no pretence to be the results of accurate laboratory experiments. More 


painstaking and less hurried practical investigation must be postponed to times of 
greater leisure. 


§1. Description of the contents of this Paper. 


Saint-Venant in his original Memoir§ of 1856 made use of the semi-inverse or 
‘mixed’ method of treatment and assumed both that the stresses perpendicular to the 
longitudinal axis of a flexed beam across planes parallel to that axis were everywhere 
zero|| and that the extension at any point was proportional to the distance of that 


* An unimportant exceptional case will be pointed out later on, § 9. 

+ We shall use ‘shear’ in the sense of a stress in which it was originally introduced into elastical 
nomenclature by George Stephenson, the corresponding strain being called a ‘slide’ by Pearson. Cf. 
Todhunter and Pearson’s //istory, Vol. 1, p. 882. 

{ Comptes Rendus, 'T. uxxxvi1, 1878, pp. 849—54 and 893—9. The paper is described in Todhunter 
and Pearson’s History, Vol. 11, Part 1, p. 191 e¢ seg. Compare also Greenhill, Mess. of Maths., Vol. 10 
(1880), p. 83. 

§ Mémoire sur la flexion des prismes, Journal de Liouville, Tome 1% (2° série), pp. 89 —189. 

|| Or, in other words, that there was no mutual stress between the longitudinal ‘fibres’ of the beam 
perpendicular to those fibres. This is the phrase used by most writers, Saint-Venant among them, but the 
suggestion of fibrous material which it conveys makes it unsuitable for use without qualification although 
it brings the physical meaning of the hypothesis well before the imagination, 
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point from the neutral plane. Clebsch* showed that the second assumption flowed 
from the first, but, as his equations are all for the isotropic case, we have briefly in 
§§ 2—4 given the required equations generalised for the aeolotropic case which we 
treat later on in the paper. 

In the second part of the paper (§§ 5—8) the general equations—isotropy 
being for the present assumed—are transformed for use with cylindrical coordinates 
and the complete analytical solution given for the case of a uniform beam with 
cross-section bounded by two concentric arcs and two radii, the plane of flexure 
being taken at right angles to the radial plane of symmetry of the beam. This 
plane of symmetry is therefore the ‘neutral’ plane as regards extension. 

Then follows (§§ 9—14) an arithmetical reduction of the analysis, especially 
in connection with the torsion due to shearing stress and its dependence on the 
angle and on the curtateness of the sector. The results, both of theory and of 
experiment, lead to investigation of the distribution of vertical shearing stress 
along the central axis of symmetry of any cross-section (§ 12), and of the nature 
of the end-fixing which is permissible with the solution obtained in accordance 
with Saint-Venant’s assumption. 

So far the discussion in connection with the particular case of the circular 
sector cross-section has been made on the assumption of isotropy (and of uni-constant 
isotropy so far as the numerical work was concerned), but in §§ 15—19 there is 
given a projective method which leads to a series of aeolotropic solutions. The 
numerical work based on this brings out the important fact that aeolotropy such 
as occurs in a wooden beam (with the grain in the longitudinal direction) causes 
a great increase in the amount of the shearing torsion. It is found that the torsion 
of the similar isotropic beam has, very roughly, to be multiplied by 24/p, where £ is 
the longitudinal stretch-modulus and w a mean value of the corresponding transverse 
slide-modulus. For very high values of H/» this rule ceases, however, to be even 
roughly true. 

In the last part of the paper will be found brief discussions of some other cases 
of cross-sections which have been partially worked out. To begin with, an investiga- 
tion is made of the algebraic solution of Saint-Venant. It is pointed out that he 
appears to have overlooked the torsional effect altogether in his discussion of the 
forms of uni-symmetrical cross-sections which arise from his solution. An attempt— 
with no great success—was made to extend his method to cases of asymmetry. 
Then follows a short description of a more general method of attacking the problem 
by aid of suitable conjugate functions. The case of the lemniscate is fully worked out. 
If there is not, as we fear there may be, some slip in the analysis, this section presents 
the interesting feature of a reversed torsional effect. Some details are given of a few 
other narnia cases. We hope to be able to give a more adequate AMEN discus- 
sion of these cases in a later paper. 

* Theorie der Elasticitat fester Kirper, 8. 74. Leipzig, 1862. 
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I. GENERAL THEORY OF THE FLEXURE PROBLEM. 


§ 2. General expressions for the displacements in an aecolotropic material 
where there are no mutual stresses between the ‘fibres’ perpendicular to those fibres. 


Let us take this longitudinal direction to be the axis of z and the axes of 
«and y to be at right angles to it. We shall denote the displacements at a, y, z 


CO ee ee 


by u, v, w and the stresses by xa, yy, 22, yz, 20, xy. Assuming uniform aeolotropy, 
and that the planes of the coordinates are planes of symmetry of structure, we 
shall use EH to denote the stretch-modulus in the direction of the axis of 2, 
fa» f, to denote the slide-moduli in the planes of «z and yz respectively, and 7, y, 
the Poisson’s ratios for contraction in the directions of « and y due to an extension 


in the direction of z*. 

The condition that the stresses perpendicular to the longitudinal axis (z) of the 
beam across planes parallel to that axis should be everywhere zero may be .written 
analytically as 


Le = yy = xy =0 


everywhere, giving immediately 


Lau 1 dv aw 
00 joy 0% 


ov 1 Ou ate 
On | oy 
The body-stress equations, also, reduce to 
0 fou . ow 
a é mF =) ah | 
av dw | 
| ee nae | es ead Oincte areas cienegemeetet vA) 
oz & if =) ; (2) 
o(E+e)+ g ees Hg phe 
Pioa\az | da) | dy \oz dy 02 


Fixing our attention on w to begin with, we deduce 


dw fi ow i Ow ‘ie Ow (3) 
02" i 0202 at oy? dz pie axoyor — Se > 


so that ie the longitudinal stretch, is linear in w, y and in 2, and may be 


written 


© = (A+ By + C’)— (a+ By + CS © (4). 


* The arrangement of the subscripts in the notation is made different from the usual arrangement, 
because so far as the aeolotropy is concerned the problem is essentially two-dimensional, 
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Integrating the first set of (1) with this value of we get 


u= —, ($4'x? + Bley + C'x) +, (4Ax’?+ Bay + Cx) 2+ F (y, 2) 
v= —9, (Alay + 3B + Cly) + 9, (Any + £By' + Cy) 2+ @ (z a)f .--(5), 
w=(A’e+ By+ C’)z—-4(Avt+ By+ C)2+H (a, y) 
where F, G', H are functions to be determined to satisfy the remaining equations 
of (1) and (2). | 
Finally we arrive at the following as our displacements, 
u= —Ty2—{k (mx — ny’) A’ + yey’ + 20} + {4 (2? — my’) A + qvy B+ eC} 2 
—4A'?+4+4A%—yyt+ Beta’ 
v= raz — {yxy A’ +4 (my — 0°) B+ ney} + {nxyA +} (my — me) B+ ny Ch 2 
—4B74+1B2-az+ye+ p’ 
w= x (x,y) — (Hm + pm — EL) lx ay" Aa Ps a Pct <(= = +) 
+ (4’0+ By+C’)z-4 icing C)2—Blat+a"yt+y’ 
 aanvaveeertaenecegteees (6), 
where 7 and a, B, y, a’, 8’, y’, a”, B” are constants of integration, and x is a function 
of « and y which satisfies the equation — 


ox exe 
Hy at Me 2 ye eM ELed 4 Cetciir cas te nesinws acess (7). 


These equations may be tested by substitution in the ‘equations of compati- 
bility.’ They are a generalisation of the equations given by Clebsch* for the 
isotropic case. The constants a”, 8”, y might of course be included under the 
arbitrary function yx, but it is clearer from the physical point of view to keep them 
separate since the linear terms of the displacements are those which must be used 
to fix the end-conditions. The terms in y give rise to a rotation about the axis of z, 
and the other notation is arranged to suit the other possible rotations which would 
be given by a=a” and B=". 

The equations (6) and (7) give the most general displacements possible in a 
material which may be supposed to consist of parallel fibres when there 1s no mutual 
stresses between the ‘fibres’ in any direction at right angles to the fibrous dvrection. 
This was the assumption first made by Saint-Venant in his treatment of the flexure 
of a uniform prism by a load applied at its free end. It is an assumption justified 
by the consideration that the stresses xz, yy), avy must all be zero at the lateral surfaces 
of the prism, and that if they are taken to be zero throughout the interior of the 
prism we can deduce a solution of the flexure problem which can be made to satisfy 
certain end-conditions of fixing and of load. It is true that in this matter of fixing 
and loading we do not presume that the distribution of stresses over the cross- 


* Théorie del Elasticité ; Saint-Venant’s annotated French translation (1883), p. 150. 
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sections at the ends is actually attained in the theory as the experience of practice 
would lead us to desire, and that we have to make use of the principle of the elastic 
equivalence of equipollent systems of loads, but Saint-Venant on physical grounds, 
and Boussinesq by analysis, have shown so conclusively the rapidity with which 
the local effects of external forces diminish with distance from the point of application 
that in any beam having a length several times greater than its breadth dimensions the 
principle may be freely applied. 

We proceed now to the particular case of the prismatic beam under flexure, and 
after that to the particular prismatic beam which has for cross-section the sector of 
a circle or annulus. 


§ 3. Particular case of the foregoing theory for the cylindrical beam under 
flexure. 

The appropriate direction for the axis of z in the case of a right prism will 
of course be that parallel to the generating line. We will take the end that is 
fixed* to be z=0, and if the length of the beam is / the free end will be z=l. 
We assume a load to act on this free end in the direction of the axis of x, and 
such that it is statically equivalent to a force W acting through the centroid of 
the end cross-section. 

Since the load acts at right angles to the longitudinal axis, the stress zz will 


be zero all over the surface z=/, that is 
(=) = (A'n + Bly +0’) —(Aa+ By +0) 1=0, 


giving A’ = Al, Bl = Bi, C=] GG 
0 


and therefore of = (Ag+ Bagge) (1-2) 0. nce. ccevepese es emeereee (8). 


Now let us consider the conditions affecting the resultant traction over any 
cross-section of the beam. These conditions are obtained by forming the equations 
of equilibrium of the portion of the beam between the cross-section in question and 
the end z=/, the only forces being the load W and the tractions over the cross- 
section. The equations are 


w= || aededy, 0=|{ pdedy, o= | { Zdedy, | 


— Wy={[(e.ge—y. %) dedy, -W(1=2)= | [v. dedy, o=|[y.zdedy | 


y being the y-coordinate of the centroid of the cross-section and the integrations 
taking place over the area of the cross-section. 


* The exact mode of fixing the end of the beam and of applying the load at the other end will be 
discussed later (§ 14). 
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— 


Consider first the integrals concerning 22. 


From 0 = | | zzdady and equation (8) we derive 
0= | {(de+ By + INGER A, 


ai o=|{(a («—2) + B(y—9) + AB + By + C) dedy, 


where (, 7) is the centroid of the cross-section, and for this equation to be satisfied 
we must have 


EGY EG a iis ene ais daw teed tmeseine kasicn (10). 
From o= | {y.2dedy, 


ee NM i i (Aare Bay Oy) dads, 
we have 
0= | [(A 2) (y-9) + B (yg) + Ae (y—9) + Ay (eB) + 2By (y-9) 
+ Axy + By + Cy) dxdy, 
and, using (10) and the property of the centroid, 
0= | {(A (=a) (y-9) + B (yg) day. 


If* therefore we take the axes of x and y to be parallel to the principal axes 
of the section through the centroid and thus make the product moment 


[[@-® Wa) aeay 
vanish, we must have 
Again from 

e w-)=|[e. ‘tadady 

and equation (8) we have 
— W(l-z) = B{ | Aer + Bry+ Cx) (l—z) dxdy, 

and using (10) and (11) 

a w-x{{A oma yt cctta: 
or A my ce eeetenes PEM ns ADSL ia ns ve Celeb oaeet conte (12), 


where 1h | | (x—x) dady, 


the second moment of the cross-section about the line through the centroid 
perpendicular to the plane of flexure. 


* It is of course not necessary to make this simplification, 
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With these simplifications the longitudinal stretch 1s now given by 
dw W 
Oey eed 
or, if we select as origin some point on that principal axis through the centroid 
which is at right angles to the plane of flexure, so that x=0, 


ow W 


eae a 
An examination of the values of u,v, w in (6) will show that such a change 
of origin will be absorbed by the arbitrary constants and the function x. 


L® ) sgl atte see Baer raat ee ee (18). 


The expressions for the displacements now become 


W J 2 F A 
w= —rye+ rE (must! — nay) (1-2) + Ble 4] yy + Bet al 


W 
v= Tea + Fa Naty (1-2) — az tye + BP <#3424), 
. W | /E£- ia NY ee Th 2 2 
w= x (x, Day | ( on He) ay? +0 (le $2) |—Bletal'y +y 
where y (, y) satisfies 
Xe 
te aot Bao Poem reeves eserecrccececcesces (15) 


To apply the remaining traction conditions of (9) we require expressions for 
the shears xz\and yz. These are 


sua Ou , Ow\ | OM av enya | (See | "1 
=m (S42) =| +X — pe neh bn y’|+B-B 


oes Bagi OU CLO pa ox W =a) “ad . 
Rom et aye ta ET me 


The first of (9) may be written 


~  (oze xy , We 
w={{(a+ (+ eo 7) 2) dedy 


oy 
from (13) and the third body-stress equation of (2), 
: /O —~ 4s) —~ 
or Ws w+{[(< («. 2) +5 («. #9) da dy 


=W+ [= (xz cos (xv) +2 cos (yv)) ds 
by Green’s theorem, v being the normal outwards, the integration being round the 
boundary of the cross-section, and in order that this may be true we must have 


22 COS (av) +2 0s (app) 30) ane amet Oss «coke (1%) 
all over the boundary. We shall refer to this as the Boundary Condition. 
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Again substituting from (16) in 


o= || gdedy 


we get o={{(z rat x — sik — —— hh) xy ~a-+a" di dy. 


The integral i | adady vanishes because of our choice of origin, and | cydady will 


vanish * <f the cross-section be symmetrical about the axis of y. This is true in 
all the cases treated in this paper, and we will make this additional assumption 
in all that follows. We thus have 


Dy A oun ” 
-{(2 ata ) dedy, 


and, transforming by Green’s theorem into an integral round the boundary of the 


cross-section, we get 


0= | Paes 1g") Vocus alder, aside eae. (18), 


an equation which we shall refer to as Condition J. 


From the last remaining of the traction conditions (9) 


—Wy= [er Yye-y x2) dxdy, 


we derive 
0 W 5 , 
e Wy= {[l> T (p, x + py Ba: (2 s = my) EP (# —3p,) wy 
W (E—7,p,-2 ; 
Hee ar Mobs) Be typ ri (a-a")e—p,(B-B")y | dady mehetaeis (19). 


This we will refer to as Condition IT. 


We have now reduced the physical problem to the mathematical problem of 
finding the function y (a, y) which will satisfy the differential equation 


and at the same time satisfy the three conditions given in (17), (18), (19) for the 
particular cross-section we have in mind. This at least is the ideal problem; in 
reality the more immediate problem is to find the particular cross-section which 
will possess as nearly as possible the special characteristics we want to investigate 
and which will at the same time permit of a suitable simple function x being 
found for it. 

In addition to our original assumption of no mutual stress between the 


* Thé integral will also vanish if the origin be taken at the centroid but soluble cases will then 


usually be symmetrical about both principal axes. 
2—2 
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longitudinal ‘fibres’ of the beam perpendicular to their direction, it is to be 
remembered that we have also assumed 

1°, that the plane of flexure is parallel to a principal plane of the beam through 
the centroidal line, and 

2°, that the beam is symmetrical about the plane through the centroidal line 
and perpendicular to the plane of flexure, an assumption which makes 1° imme- 
diately satisfied. 


§ 4. The case of isotropy. 

It will be convenient at this stage to postpone the aeolotropic discussion and 
to summarise the foregoing for the isotropic case where y,=7,=7 and p,=~.=p. 

We have taken the longitudinal axis of the beam to be the a of z and the 
plane yz to be a plane of symmetry, the flexion taking place in the plane za. The 
fixed end of the beam is at z=0 and the loaded end at z=l. 

The displacements* are 


Ww 
w= —Ty2 +5199 (2 —y') (1-2) +312 — Fe] —yy + Beta 
parent ayy nay (1—2) + ye 02+ B Wy ps't (20), 


w=X (x, y) an [wv (z—12*)+ay?|—- B’ae+a"y+y’ 


OX eae 
Oa Oy?” 


and 7, a, B, y, a’, B’, y’, a’, B” are arbitrary constants. 


where 


In addition we have the following conditions to be satisfied. 


(i) The Boundary Condition : : 
2 COS (xv) +42 BOS (277) 0 ott cc scan ee sees Mae (22), 
h We ox W Lae at 2 " 
Wee a=) —TY +0 — a7 (Se +(1—3n)¥')+B-B 
~ oy W ‘ 
and gen p {eX — 7247) cy ata 


(ii) Condition I: 
0= i (x(a Lal) Ween. nyse (23), 


the integral being round the boundary of any cross-section, and 


(iii) Condition IT: 
W : 
— Wy= | {[rc (a? + y") + (= X— y X) — Ry 2am ey 


W 
+7 (1 —iy) y’—(B—-B") y | Boe eos. veo (24). 


* Compare Love’s Elasticity, 2nd Edition, Chapter xv, Equations (12), p. 319. The axes of coordi- 
nates and notation used in the present paper make comparison easy with the results given in Chapter xv. 
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II. Cast,or A BEAM HAVING AS CROSS-SECTION 
A CIRCULAR OR ANNULAR SECTOR. 


§5. The equations suitable to the case of the circular sector cross-section with — 
isotropy. 


For our problem of the sector-shaped cross-section we must now transform 
our equations to cylindrical coordinates 
z,7, 6. Let O, the centre of the circular 
ares be our origin, and let the axis of y 
be the line 6=0. We shall take the 
sector to be bounded by the ares r=a@ 


and r=a,, and the lines 0= + 4y. Ota 
The boundary condition 1 
“— — ! 
az cos (xv) + yz cos (yv) = 0 RY, 
: : i ~ ma X 
is, on substitution of az and yz, 


0 0 
sa cos (av) + a cos (yv) 


=+ (y 00s (av) — 2 cos (y»)) + jy [C92 + (1 — 41) y") €08 (ev) + (2+) ay cos (y>)] 


— (8B —B”) cos (av) +(a— a”) cos (yv)... 2.006000 (25), 

and we can now obtain the conditions that are to be satisfied on each of the 
edges of the sector. 

1°. Over the external arc, where r=a and —ly<6<+4y, the boundary 
condition becomes 
W 
EI 
+(2+7) a’ sin 6 cos’ @]—(B — B”) sin 6 + (a — a”) cos 6, 


(*) =7 (a cos @ sin 6— asin @ cos 0) += [(Sna’ sin’ @ + (1 — 4) a cos’ A) sin 6 
Tr=a 


() mn (?+4y) sin 0 +2 sin 36]—(8—f”) sin 0+ (a—a”) cos 6...... (26). 
2°. Over the internal are, it is easily found that 
(| z es (2+) sin 0+2sin 36]—(8—”) sin 0+ (a—a”) cos 0......(26)’. 


3°. Over the radial edge, a,<r<a, and 0=4y, the condition is, in the 


same way, 


a Wr’ 
("2B ecyy =? + gy [t= ¥0) 2088 +4 008 87] 
—(B— 8") cosdy —(a—a”) sin Fy...... (27), 
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and, 4°, over the other radial edge 0= —4y, 


won = — [(4-4 5k 3 3 
Boe, Ue EI [(4—4%) cosdy +3 cos By] 


—(B- 6") cos 4y + (a—a”) sin dy...... (27), 
the only difference between 3° and 4° being the sign of the term containing a—a”. 
But if y is to be one-valued, the conditions 3° and 4° cannot both be satisfied unless 
a=a"’. It will be evident after we have found y that the physical interpretation 
of this is that, in a beam which is symmetrical about the plane of yz, there is no 
lateral slide when a flexing load is applied in the principal plane za. 


§ 6. The solution for x sutable to the case of the sector. 
We shall now assume*, as our solution of 


OX Oe 
0a? tay =), 


Vas (4, em 4 =) sin m0 +a,r sin 0+ a,7° sin 20+a,7° sin 30 ...(28), 


an 
mm 


and proceed to determine a,, @,, @,, A,, A,,..., B,, B,, ... and the appropriate values 
of m from the conditions obtained above. ; 


From the boundary condition over the straight edges 


f) B 
(% => Sah ee =) mcosémy +a, costy+a,r.2 cosy +a,7". 3 cos 3y 
00) 9=44 


<a ee sey fee 3n) cos y+ 2cos8y|]—(B—f”)coshy ...... (29), 
which will be satisfied only if 
my _ 1 
5 = (24 Bais: 


where 7 is zero or a positive integer, 
T 1W 
i 
and a,= —(B—f’”), 


d= 
Hence we have 


~ 2cosy’ 


on Be Q VIN Och aco ae wi A: 
x= = (aie + wat) sin m0 — (8— 6 )a sin +5 meh sin 20 


* We have to make x an odd function of @ to satisfy the equations of §5. Physically this means that 
the longitudinal displacements in the beam are symmetrical, with change of sign, about the axis of 
symmetry of the cross-section. 
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Applying now the boundary condition over the arcs, we have 


(4) = % (Angmar —m pt =) sin m0 
Or /r=a 4=0 a 


W cos 4 : 
—(B—B") sin CS eenes asin 20+ ar E +(4-—4”) =| v sin 30 


ene | +4n) sin 0+2 sin 30 | — (B—B") sin 8, 


or 


m+1 


a2 B 
> (A,,a" 1 —-— |) msinm# 
i=0 a7 


ids a TO Wa? cosdy . 
= = ik = of 1 1 oY ‘ 
7 (at) sin @ Segal 20 EI (4-47) a 30 ...(31), 


along with a similar equation in a, for the internal surface. 


These we may write for brevity as 


0 


a 


= C,, sin mO=Se, sin pO — [p= 1, 2, 8] ersreeeee. (32), 
: ; 
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where c, is known and C,, contains the quantities A,, and B,, which are to be found 


We can now determine the constants C,, by Fourier'’s method and thereafter 
obtain A,, and B,,. Multiply both sides of the equation (32) by sin m@ and integrate 


between the limits —3y and }y. Then 


2 2 
Cn | "sin? mOd0 = ze, | ” sin pO sin mOdé 
~4Y 


= 32% ie 2 {cos (p —m) 0 —cos(p+m) 6} dé, 


and trae She pinkde=nir em) 
Y p oa ee parm 
3 4p cos dp . 
=3 ceri aren Newer tater eeici’ t, (33). 


Replacing the notation C,, and ¢,, we have, after a little reduction, 


B,  W 3+2n coss4y 87 (—1)* 


MA Neer qnt! “Ei nm a) y = 1) ee (m esa ah 


may lL 25 oY (= 


1) a 


EL on — 
Ae. Boned VW. i, cosh a 87(— 1) 
MA mM Ne ht = AL m — (- l)a, 0 ~ (m?—4)- y Uy 


_3W1—2n costy ae 
ay ARRAN 
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the solution of which is 


oie ts 18 isin) Mace 


A,,=-— 
zi (m?—4)my a@™—a,™ 
+ (Lay pey Ceeeae (Sen) ee 
KI my Hoe Al —9 qe — iP a 
B or ( dt by wad ues (ae oe ats ). 
m (m* — 4) my qm — on. 
i ( a 1) W cos sy 3+ 2n uy 3 e Mt 2m) ots We (an a5") 
KI my Ui fie AL <2XG a — oo 


We thus have y and therefore the displacements wu, v, w expressed in terms of 7 
I. 


and the arbitrary constants a, B, B”, y, a’, B’, y 
W 
w= —ty2+ Fr[S (ne — ny’) (1-2) + $12 —42)]—yy + Beta! 


W 
v= Tae + Fay gay (L—2z) — az + ya + B’ 


pth : & S Sr ( 5 By ant =n Go. id = 81 ( et 1) iit Asie (a pets da i “| 
red (m a 4) my qm — Ge mae (m Ce: 4) my ae” — ENohe ve 
= W costy (8+2y 3(1—2m)\ am? - an on 
A (= Dy EI my (= — | m Nan 9 am a ri gs 
Say W cossy /3+2y 3 (1—2m)\ a™*ay™"* (a™*—a"*) 1 | sin m0 | 
a6 EI my \m?-1 = m-9 a” — a," co | 
W cos 
SLe faye i ae cee 27 
(B B") rsin 0+ 5 een *sin 20+ spr | 2+ (4 A) Sosa |"? sin 80 | 
ee 
— ay le (2— $2) + ay"] — Bla t aly +7 | 
PDE ERRTIA HAO ARS 5. (35), 
where x=rsin 0, y=rcos 0 
and m=(2i+1)". 


We note that the displacement in the direction of the flexing force, u, is 
symmetrical both in sign and magnitude about the radius of symmetry, while v 
and w have both alternating symmetry. 


§ 7. Application of the Conditions I and II and deriwation of the value of 
the torsion t. 


We have now to apply the Conditions I and II, equations (23) and (24) of § 4, 
and will thus determine two relations between t and B—£” from which these may 
be obtained. 
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(1) Condition I, that is 


[x COSY E) S30 Ayana canon so; se sucepawe ests (36), 
the integration being round the contour of any cross-section. 
Taking the four parts AB, BC, CD, DA (see Fig. 1) separately, we have: 
1. Along AB, cos(yv) = —sindy, and 


B a 
[_ xe2s (yv) ds = —sin 4y [’ Xpaty OT 


2°. Along CD, in the same way, 


D ‘ag B 
| x cos (yv) ds = sin by | Xpa-y (ar) = -| x cos (yv) ds. 
6 a 4 ee A 
3°. Along BC, cos (yv)=cos 8, but x is odd about 9=0, and therefore the 
integral vanishes. 
4°, Along DA, the integral again vanishes. 
Thus Condition I is identically satisfied. 


It is easy to see from (36) that Condition I is identically satisfied whenever the 
section has uniaxial symmetry about an axis perpendicular to the plane of loading. 
For in this case x changes sign as we cross the axis, and, if P and P’ be corresponding 
points on either side of the axis, an element ds at P will correspond to an element ds 
at P’, cos (vy) will be the same at both points, but y will change sign as containing 
only odd powers of x, and thus the two elements of the integral in (36) just cancel. 


(II) Condition II, that is 


es ene Ore, ox 
“tives 99) 


— Fay (2+ bn) y+ a py tan ¥—(B-B") y | dey. 


When the integrations are carried out, we obtain 


. [lex- A) dady = — pr 4 (a —a,') + (B—B") » 3 (a —a,") sin $y 


4a —a,’ sin OW ies al aR : 
= ae OY Foe (ata) —40) sim y+ 48in By}--(37), 


and, transforming the surface a by Green's theorem, 


[[(e2 eA ay X) dedy = [x (x cos (yv) — y cos (av)) ds. 
Taking the line ee in four parts, as in the previous case, we have : 


1°. Along AB,  cos(yv)=—sindy, cos(#v)=cos$y, and 


B a a 
[, x (« cos (yv) —y cos (xv)) ds = -| Xonay (7 Sim’ by +7 cos? sy) dr= -|' Xgapy 7" 


3 
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2°. Along CD, in the same way, 

[ox (a cos (yv) —y cos (xv)) ds = + Xonay? (— 41), 
and therefore the sum of the values along AB ah CD 
— -2[" Xpngy TO 
3°. Along BC, 
[x (x cos (yr) — y cos (av)) ds =|x (asin 0 cos @—acos 0 sin 8) ds = 0, 

and, 4°, along DA the integral again vanishes. 

Thus the whole integral, P, is simply —2 i oes rdr. 


Substituting y from (30) we find : 
Cae —* pat ‘ fie | 
°—— sin 4my 


m 


P= = 9 EE As (a? — a"*?) + 5 B 


Nezo m+ 2 j-0m— 2 gm? C, 
—(B—B")4(a—a,) sin by +74 tan y (a —a,') 
1 Woe cos 4 ; ‘ : 
a 15 AL (3 +(@—37) a gaye —ai)| --.(38), 


A,, and B,, being given by the formulae of equation (34), and on the substitution of 
this value of P in (37), of the values of A,, and B,, from (34) and the value of J we 
obtain, remembering that # = 2u(1+7), after some reductiotis the following equation 


for T: 


sa) 64 (a™? _@ pas 
cay iB gcas 4 t a? > 0 
| (a*—a,') (tan y—y) + 0 (m—4) (m+2)my a ae” | 
ef a) 64 a me: Ce | 
7=0 (m’ ae, 4) (m — 2) nvy Ores =a es 
~ BL Zo(m+2)my\ m1 m —9 ioe 
= 3 oy ee Ce - 
o(m—2)my\ m1 m?—9 a” — a" 


+ #5 (1 —2n) cos $y tan By (a’ — a,.’) 


-Hs 2.cos hy te Are CAG: anh Ds 9 


sin Ct — Ay’) (a — Ay’) _ 
4 ( +n) ( i sin sy ree 2(1—2n) sin dy (a@—-a, 
The value of B—£” still outstanding will flow from the nature of the ae 
conditions. 
As we shall see in § 9, it is this quantity 7 which gives the torsion in the beam 
due to the asymmetrical distribution of the shears in any cross-section. 
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§ 8. The solution in the case of the complete or non-curtate circular sector. 


For the case of the complete sector, that is when a, is taken to be zero, the 
solution of §§ 6 and 7 simplifies into the following : 


T 


x=24,,7" sin m0 —(B —B”)rsin 0+ aoe rv’ sin 20 
W cos 4 
“ ig ag ZV loses 
| 3 KT [a+ in) ale sin 30 eeeeee (40), 


where A,,=— : AS 
m—1 mM —9 


87(-1)' 1 ,(-l)! W costy 38+2y 3(1-—2n)\ 1 
™ (mm? —4) my a EI my 7 


and m=(2%+1)", %7=0,1,2,... 
y 


and the torsion is given by 
64 2 i! 
oy ‘ 
r| Se a y aa (m? — 4) (m+ 2) =| 
ue: 2 cos by (Het e0 _ eC 


~ ET Lio (m+2)my\ mi—1 a —9 


> 


00} 
Fn 
il 
+ 

3 
~~" 
7 Gea 
fad 
| 
|e. 
215 
fess 
Sees 
=} 
bj 
SS 
ae 


++ (1 —2n) costy tan gy —2 (1 — 2n) sin $y —¢ 


Ill. Numerican WorRK IN CONNECTION WITH THE TORSION DUE TO FLEXURE. 
Discuss1on OF SHEAR DISTRIBUTION AND OF TERMINAL CONDITIONS. 


§9. Torsion of beam with flexure. 


As stated at the beginning of the paper the torsion of the beam about its axis 
is due (i) to the resultant torsional couple which arises from the shears in directions at 
right angles to the axis and (11) to the anticlastic curvature caused by the inequality 
of the longitudinal tension and consequent inequality of the Poisson contraction over 
any cross-section. For each of these effects to cause torsion it is of course essential 
that the cross-section be asymmetrical. 

To determine the amount of torsion let us take the vertical displacements of 
points on the horizontal axis of symmetry at the fixed end and at the other end 
of the beam. These are given by the equation (20) 


w= — ryt gy lh 2) (wy) FMP Be yy t Beal 


Firstly, at the fixed end of the beam the displacement, (i) at the outermost point 
of the axis of symmetry, namely (0, a, 0), is 


W : 
Ug = ry Lyla —yat+a, 
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and (ii) at the inner end of the axis (0, a, 0), 
W 
tL Rise EI $yla, —Y% ae oe 


giving Ua — Ua, = — se dyl (a —a,’) —y (a—4,). 


The first term arises through the anticlastic curvature and the second, being 
merely a rotational term affecting the whole beam, depends on the particular method 
of fixing the end of the beam. When we consider the torsion of one end relative 
to the other, it will cancel out. 

At the end where the load is applied the corresponding displacements are 


(i) at (0, a, 2) Ud = —Tal+ Ty abaya Bla 


and (11) at (0, a%, /) Ua, = — TA,l + ue £P—ya,+Bl+qa’, 


giving Ual —Ual = —T1(a—a) —y (4-4), 
the terms El 41'+Bl+a’ disappearing since they are really the expression for the 
flexure. 
: x. Ug tee a 
Taking p= i tyl(ata,)—¥y 
Ua! — Ubay, 
and wi = ae = ql rvs 


as approximately giving the small angles of rotation at re 
k =a 

the two ends of the beam, we have the total torsion - — * 

of the loaded end relative to the fixed end about the 


longitudinal axis of the beam 


I 
i 


Ww 
=, -— Y= — th + ay 4yl (ata) SaSeesod (42), 
7 Direction of 
the direction of the torsion ,—~, being as shown in Flexure 
the diagram of the free end of the beam (Fig. 2). Fie. 2. 


In connection with the anticlastic curvature it is to be noted that it is zero 
at the free end, there being no stretch there, and that, accordingly, there is no 
distortion of the surface at the free end except for the Saint-Venant distortion in 
the direction of the beam’s longitudinal axis. Nevertheless if we were to clamp 


* At first sight this expression for y looks strange in that it would appear that the angle of torsion 
due to the anticlastic curvature would be infinitely large for the limiting case of a rectangular beam 
which is given by taking the curtate sectors to have infinite radii a and a, but such that a —a, remains 


finite. In taking the limit of #1 (a+ a), however, it is to be remembered that we must take the limiting 
expression for Z as well, and J being } (a‘— a,*) (y — sin y) it is seen that the angle of torsion is proportional 


to =o —, which gives zero angle of torsion for the rectangular beam as it should do. 
(a? + Gy") (@ + a) 
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the fixed end in such a way that the ends of the axis of symmetry (0, a, 0) and 
(0, a, 0) were fixed, the anticlastic curvature of the beam at that end would give 


rise to the torsion aq (a+a,) at the free end. If the fixed end be clamped in 


another way, the visible torsion due to the anticlastic curvature would be less and, 
it is clear, might even be zero. This would happen, for example, if the clamping 
were such that it fixed that element of the axis of symmetry, «=0, which, after 
the curving due to the flexure, remained in a direction at right angles to the plane 
of flexure. In the practical case of the beam which approximates both in length* 
and cross-section to the shape of the section of a propeller-blade, it is found that the 
anticlastic curvature gives rise to a rotation about half that due to the shear torsion 
but opposite in sign. 


It must be borne in mind that the anticlastic torsion depends on the nature of 
the ‘fixing’ and the methods of clamping used in practice differ very widely from 
the theoretical means of fixing provided by the constants a, B, y, a’, 8’, y’, a”, B” of 
our displacement equations. These constants occur only in the linear parts of the 


equations so that the most we can do with them is to fix certain points and certain 


; a(t : : 
coefficients of the scheme ae ‘ 3 at points on the plane z=0. Remembering that 
we have found that a=a”, and that there is an equation for (8 — 8”), we see how easy 


it is to exhaust the stock of arbitrary constants. Thus, if we try to fix an element 


: ; Ow ow ou/ ow 
at the point (x, y, 0), we find that making w=v=w=0 and ae Wy eae (- | =='() 


gives all the arbitrary constants definite values. Similar difficulties arise with the 


method of loading, where the assumptions of the theory are even more difficult to 
approximate to in an actual experiment. 


§ 10. The torsion of the non-curtate sector for various values of the angle vy. 
1%. y very small. 
To get an approximate value of + when y is small we can expand the trigono- 
metrical functions and the three series of (41) in powers of y. 
zy 1 


If n= iat (m aa 2) ny (m’ i n’) 


f= De2s 3], 


then S,= > 


p=0 tae i 
% ae {at {Us ¢ bisa Bon (2) } j 


WhETC Ue p1= 2 nllyis nUopsg = nly + 2” *, and .u,=1. Further 7,= - (Qi 1)” and its 
values are known and tabled, while we ourselves have tabled ,u, from qg=4 to 31. 


The series, however, does not converge rapidly and can only be used for y small. 


ao 


* For the torsion is proportional to the length. 
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In the latter case, omitting much troublesome algebra, we find : 
T= He [de + 12 + (084,033 —*168,066m) y — (254,632 + 032,403n) ¥" 
+ (027,637 —*198,882n) y’ — (-084,951 —130,7767) y‘ + ete. ]...(43 a) 
- ih [1 +°042,017y — 262,733 — 022,083y' — 052,258y'—ete.] ...(43 b), 
when we substitute the uniconstant value 4 for Poisson’s ratio 7. 
The limiting value as y—0 is therefore 


Wa 
cmd x 28333. 


; : W 
2°. y=3'. The approximate formula gives 7= 4 “ %-28481, the accurate value 


El 
has the coefficient *28504. 
X 28482. 


3°. y=6. Using the approximate formula (43 b) we get r= it T 


The accurate value of the coefficient is °28577. 


: : : Wa 
4°. y=12°. The approximate formula (43)) gives 7 = Fy * 28031, the accurate 


value being *28461. Clearly, except for very rough work, the approximate formula 
cannot be used beyond about 5° or 6°. : 
In aconsiderable number of cases the series were summed and the results are stated 
below. This number could easily be extended, but we found the labour greater than 
calculating by modern processes the series to even 80 or more terms. Great care has 
to be taken that the right limits are found for series and terms which give infinite 
values. These infinities always cancel, but the difference of the two infinities may 
provide a finite term. Interesting summations bearing on the theory of numbers 
constantly occur. 
5°. y=45°. This case is analytically of some interest for 
64 
i=o (M+ 2) my (m (m ? 


“M8 


=4) = < (Slog, 2-ta+ pce 


the ten-figure decimal being provided by Glaisher’s sum of 3 : see Proc. 


o (40+ 3)? : 3)" 
Lond. Math. Soc., Vol. vim, p. 203. 
6°.. y=90°—a beam with right-angled sector as cross-section. 
In this case m=2(27+1), and the series are all expressible in finite terms. 
Thus, omitting all the laborious algebra, 


‘ 64 

=—2 Adisea) 9 
nh (1 + 2) my (m? — 4) —tany+y=37 = log, 2, 
Pan (m+ 2)my (m? as 9) = ee (6 log, My —7), 


: 1 “hoa 
i=o (M+2)my(m?—1) 1807 
Substituting these in (41), we obtain 
v Warne (112 log, 2+ 4 — 267) + (16 log, 2+ 52—187) » 
EI 107° — 120 log, 2 


(16 — 32 — 22 log, 2). 
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7°. y=180°—a semi-circular solid cylinder with flexing force parallel to the 

diameter of the cross-section. By similar summations we find: 
_ Wa Sy hia 2+ 6y 
"EL w-8 15 

8°. y=270°—a solid circular cylinder with a quadrant cut out of it, the flexing 

force being perpendicular to the bisector of this quadrant : 
Ps J2 _ (4556 42352 log, 2 — 5467) + (3068 — 336 log, 2— 3787) y 
Leet eS 15 20° —7 —8 log, 2 

9. y=360° 

central axis, the flexing force being perpendicular to this plane cut. 
Lee 8d (8+ 9n) 

TEI ” 25 (90° — 64)’ 


Collecting all our results we have the following table : 


a complete solid circular cylinder with a plane radial cut up to its 


We find that 


Tasie Ll. Torsion in a bent beam which has for cross-section a non-curtate sector. 


n=} n=% n=4 
f Wa = - - i 
ea Mig | EL ee Wa Total torsion Wa_ | Total torsion Wa Total torsion 
Diels Wael. / Hi |. Walk rf EI WajEl 

0° 066,667 + °866,667y | 283,333 | — 158,333 | +355,556 | —-188,889 | -500,000 | — 250,000 
Di 070,565+ °857,911y| °285,043 | —-160,043 | 356,535 | —-189,868 | -499,520) — :249,520 

4°-5476 | Interpolated max. y = 4 — — 356,652 | — 189,985 = = 
6° 073,458 + +849,257y | °285,772| —-160,772 | -356,544 | — 189,877 | :498,086 | — -248,086 

6°-9426 | Interpolated max. y=}| °*285,816 | — 160,816 — — = =o 
Ip 076,642 + -831,8607 | °284,607 | —-159,607 | °353,929 | — -187,262 | -492,572 | — -242,572 
ES" 076,772 + °814,098y | :280,297 | —-155,297 | 348,138 | —-181,471 | -483,821 | — -233,821 
24° 074,327 + -°795,822y | °273,283 | — 145,283 | -339,601 | —°172,934 | °472,238 | — 292, 238 
30° 069,762+ -776,990y | °264,010| —:139,010 | +328,759 | —-162,092 | -458,257 | — 208, 257 
45° 051,839+ °728,126y | 233,870 | —-108,870 | -294,548 | —-127,881 | -415,902 | — 165,902 
54° 038,645 + -698,488y | +213,267| —-088,267 | :271,475| —-104,808 | -387,889 | — 137,889 
60° 029,643 -679,077y | 199,412 | —-074,412 | -256,002 | — 089,335 | -369,181 | —-119,181 
75° 009,098 + -633,707 | °167,524 | — 042,524 | -220,333 | — 053,666 | 325,951 | — 075,951 
90° —004,459 + pe eg 144,580 | -- -019,580 | +194,259 | —:027,592 | -:293,618 | — 043,618 

104°-8863 Interpolated min. y=4)| 136,646 | — 011,646 — — — — 

106°-1906 ra - n= = — 184,139 | — 017,472 — — 
108°:3252 n=4 — —_ — — ‘278,781 | — 028,781 
120° 006,222 + "558 56557 | *145,886 | — -020,886 | -192,440| — -025,773 | +285,549 | — -035,549 
150° 082,335 + 584,380n ‘228,430 | — :103,430 | -277,128| — 100,461 | -374,525 | — -124,525 
180° 224,047 + *672,141y| °392,082 | — -267,082 | 448,094 | — -281,427 | -560,117 | — -310,117 
210° ‘407,389 + -798,737n | -607,073 | — °482,073 | -673,635 | — 506,968 | -806,758 | — 556,758 
240° ‘597,136 + -929,6507 | *829,549 | — °704,549 | -907,019 | — :740,352 | 1:061,961 | —-811, "961 
270° 755,377 +1:028,203y | 1-012,428 | — °887,428 | 1-098,111 | — 831,444 | 1:269,478 |-1 019, 478 
300° *852,715 +1:066 mee 1:119,321 | — 994,321 | 1:208,190 |-1 O41, 423 | 1°385,928 |-1- 135, 928 

| 319°-4474 Interpolated max. 7 =} | 1:138,590 -1-013,590 -- — —- 

318°:2777 Fp x n=4 — — 1:226,331 ~1-059, 664 a = 
316°°3202 i .f n=4 — a | = 1:402,199 |—1:152,199 
330° 875,056 +1-031,731y | 1-:132,989 |-1:007,989 | 1-218,967 |-1 052, ,800 | 1:390,922 |—1:140,922 
360° 824,927 + +928,043y | 1:056,938 | — -931,938 | 1-:134,275 | — 967, 608 | 1:288,948 |-1:038,948 


The total torsion =anticlastic torsion minus shearing torsion, and it tends to move the angle edge of 
prism in the opposite direction to the flexing force. Since the total torsion is always negative the angle 
edge of prism will always move in direction of flexing force. 


Wa 


Value of + / 
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The curve in Fig. 8 shows how the nature of the torsion changes with y. It 
provides a kind of physical measure of the asymmetry of the circular sector, for, 
as we have seen, the effect is wholly due to the distribution of shearing forces over 
an asymmetrical section. Remembering the direction of the torsion as given by 
equation (42) and Fig. 2, we see that in the case of a wedge-shaped sector of small 
angle, such as the blade of a sword held horizontally*, the torsion will be such that 
the back twists up with a downward flexing force. A first maximum is reached 


1:40 


EI 
° 
3 


2 
a 
ro) 


Base of Total Torsion =z 


oi Base of Total Torsion n=% 


0:00 
or 20° 40° 60° — 80" 100° 120° 140° 160° 180° 200° 220° 240° 260° 280° 300° 320° 340° 360° 
90° 110° 130° 150° 170° 190° 210° 230° 250° 270° 290° 310° 330° 350° 
Angle of sector (y) 


Fic. 3. Torsion due to Shearing Forces in a flexed beam the cross-section of 
which is a non-curtate circular sector. 


between 4° and 7° according to the value of Poisson’s ratio. As we increase the 
angle of the sector the torsion diminishes in magnitude until with y equal to about 
105° the value of r is a minimum. From consideration of asymmetry this is what we 
might expect, for with y=105° we have from various points of view the sector which 
most nearly approaches symmetry. With angles larger than this angle of approximate 
symmetry the torsion of the beam increases again. With the increasing departure 


* A rough experiment was made with a carving-knife to which a long light pointer was attached 
to show the torsion when a weight was hung on the end of the blade. The torsion was of course 
extremely small, but was seen to be in the direction given by the theory. 
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from symmetry as y becomes larger, the amount of torsion continues to increase until 
a second maximum is reached about 316° to 319°, from which value it falls slightly 
until y= 2a which is the physical limit to the application of the formula. 


§ 11. Numerical examples of the non-curtate sector. 


The total torsion in a beam of length / is given by 
V 
T= te x the constant tabled in Table I 
V 3 
= sate o x the tabled constant 


3 
= 0 x > x the tabled constant, 


if 3 be the deflection of the centroid of the end section. This is on the assumption 
that the length / is so great that we can neglect the linear term s,/ (see equation (50) 
mfra) which comes into the formula for the deflection owing to the shearing forces, 
an assumption which we shall investigate later in the paper*. 


TasLe II. The deflection and torsion of steel beams of circular sector 
cross-section with a load W lbs., for varying size of angle. 
Frencil of beam = 20 inches, breadth = 1 inch, stretch-modulus = 30,000,000 Ibs. per square inch. 


Shear, torsion per 
bicegid Deflection sav inches | Shear torsion=7l degrees Sal ae ase 

ily W x 802-53 W x 97°9356 0°-122 
3° 29s x 3°6416 0°-122 
6° a Sa TT x 0°45646 0°:123 
12° x 0465 x 0:05686 0°-122 
18° x _0°138 x 0:01662 0°-120 
24° x  0:0586 x 0:00688 5 O17 
30° x 0:0301 x 0°003,415 0°:113 
45° x  0:00908 x 0:000,9123 0°-100 
60° x 0:00393 x 0:000,3368 0°:086 
90° x  0:001,246 x 0:0000,7741 0°:062 
120° x  0:000,5789 x 0:0000,3629 0°-063 
150° x 0:000,3357 x 0:0000,3295 0°:098 
180° x  0:000,2264 x 0:0000,3815 0°:168 
210° x 0:000,1707 x 0-0000,4453 0°:261 
240° x  0:000,1407 x 0°0000,5016 0°:356 
270° x  0°000,1245 x 0:0000,5417 0°:435 
300° x  0°000,1165 x 0°0000,5604 0°-481 
330° x  0:000,1136 x 0:0000,5531 0°:487 
360° x  0:000,1132 x 0-0000,5141 0°-454 


Norn. <A positive torsion signifies a torsion in the opposite direction to the arrow in Fig. 2, :p. 20. 
In Table II we have as a particular case the data for a steel beam of length 
20 inches and radius 1 inch, the load applied being W. Taking # to be 30,000,000 lbs. 


per square inch and calculating J from the formula a (y —sin y) we get the deflection 


* See § 14. 
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and the torsion of the beam. For another beam of length Jl’ inches, of radius 
@ inches, and stretch-modulus E’, the deflection is got by multiplying the deflec- 


tions of Table II by 3,750 and the total torsion by multiplying the tabled values 


L? 
Ea {4 


Ig 
by 1,500,000 are sa 
§ 12. Distribution of vertical shear along the central radius of a cross-section. 
To get some idea of the distribution of shears which gives rise to the torsion, we 
will consider the distribution of the vertical shear along the central radius of any 
cross-section, that is the shear zx along the line #=0 in any cross-section*. We shall 
confine ourselves to the case of the complete sector. 


From equations (16) and (40), 
a oy W : : [ 
c= p| ry +X — 7 (ne + (1 dn) +88" 
= {= TY — ate +(1—4n) y’)+B-f’ 


@ | , 
a E A,r” sin m0 — (B — B”) r sin 8 ort r sin 20 


W 
ert +(b-39) Sad) sin 30 || wits | ieee eee (44). 


Along x=0, or 6=0 we have 


aoe TY W ¢ cos 4 2 
20 = H|-- AAs n) YP +LApmy +t + ard G-ay) ) y 


if 1 W ; cosdy\ 
=H \7Y (say + ay(-$+3 8 cos 3 ay 


Sy eat 1 yeu! Wee Ha Adee 33 8 1 jy} 
ai (m*—4)y a" KI y \2(m-1) 2(m'-9)) a” J ‘ 


when we put »= % and insert the expression for A,,,, 


We [fh pho fal 1 /cos dy y 
RT. (7 5 Te (cay) +gleeee—1)% 


/ Wa 
i 


Die} 
fo) YONG 7. 3 ) Pres ‘ 
po MepeS ee Sa a, Dee ma es a aie 
+3 araa tom by (5 Geena) — aterezay) ya) tae ete 


where m=(27+ 1) a 
Y 


* Shears other than in the vertical direction will of course have appreciable and in certain cases may 
be of great effect on the value of 7, but the labour of calculation would be considerable except in the case 
we have taken, and that, after all, is the one likely to be most instructive. An important result also flows 
from this central radius shear: it enables us to determine (see p. 44) the shear at the centroid which 
provides the linear term of the deflection with Saint-Venant’s form of the terminal fixing. 
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The form of equation (45) allows us immediately to compare the distribution 


of shear in beams of various angles, but of constant length and radius, when the 
CLARE nee : 
deflection = =,, is kept constant. This is the suitable method of comparison, and 


BBL 
in what follows we will only consider zat / (1 DD , that is the expression in the double 


: : ; the vertical shear on central axis 
brackets, it b rtional to — ; : 
rackets, it being proportional to qcaanten 


It is useless to attempt to see the nature of the distribution of shear from the 
equation itself, and we have immediately to resort to numerical calculation. 


Wa? 
in i) for every tenth of the 
radius of the sector and for sectors whose angles are 5°, 12°, 30°, 45°, 60°, 72°, 90°, 
120, 150°, 180°, 270°, 360% 


In Table III below are given the values of “a / ( 


For all these values of y the series involved in equation (45) may each be summed 
in finite terms. For let us consider the series 


ater 


i-0 mM —/p’ \a 


a series which for p=1, 2, 3 gives the three series concerned. 


if y m-1 | ae) i} 1 y m~ti1 © ‘eee 1 LYN By | 
ae (") eae Op E m—p (¥) baad) ~m+p (a) a) 


BPE d\. 3 DO Nee i a 1 ne *" 7 
2p (4) = cael ey AG Lele cad) 


y ¥ 
oy y b ae | a y m—p-1 ue f z) ¥. ¢ i a [; (“yr a j Al 
bP | (4) > [; (2) rg (J t) Sa hO NG ee ie 
Ban! (2i41)"-p-1 
Bay SON eet eta) (ay u) 
2p ( , {; Ee (7) koi ( 


z , (+1) = +p-1 


ea ag) 


the interchange of the integral and summation signs being clearly justified, 


y 


: 2 ae 
ioe oc tO 
apa y a Rey As u\ - ae a | ( 
2p G Te € (z ae T(E) Jere (46), 
eh op 
e a 


: : : : 2a is : 
and both of these integrands being rational functions of af z , when - is an integer 


(as in all the cases considered), we must be able to express the series in finite terms 


by evaluating these integrals. In general they are incomplete B-functions. 
42, 
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For the sectors whose angles are 5°, 12°, 30°, 45°, 60°, 72° and 120°, it is easier to 
use the series as they stand than to find the sums in finite terms and calculate from 
them. For these cases the equation (45) becomes : 


—~ 


for 5°, Zi f(u-pr)= + 0010917 +: ‘000959 (%) — ~ 002668 (2) + 00028 (2) ae 


for 12°, <n] ar)= a: 0063584 4: ‘oos713 (4) — 007971 (4) 


104 


+-o00701 (4) ~ -000247 (4) 4000128 (4) ep 
O Ob O 
Phat haae. Ah y\" y\" 
for 30°, za] bpp )= + 7040842" +-045753 (4%) —-044065 (4 
17 29. 
+ 001400 (¥) ~ -000426 (4) Bes 
a a 


for 45°, Zafar )= + 0968727 4° ‘176777 (4) ~*176109 (4) 


27 


11 19 
+ 001295 (4) — 000210 () + 000073 (“) ee 
Ob Oo a 


for 60°, zi |(w py )= + 1994124 —-ossar (4) 4760556 (2 2) logy! + ‘ooosse (4) 


\ 14 


/ ; 20 26 32 
+ 000286 (2) = -oo0221 (4) + 000152 (¥) ~+000108 (4) tea. 
O a O a 
5 “ a. Wa . y 5 oc y . 
for 72°, <n |(u zr)= + 38699 %—-452254 (“) 
; 6 Pash 16 
+ [4 (2807814 + 000780 (“) ++(000706 () ~ -ooo4si (2) +e) 


for 120°, ca] a) = ~ 437658 4 —-187500 (2) 


} 
+ 000209 & i 
) 


~-000134 ( a + 000105 Z) + 000070 (2 


“ = 000050 (4 © + 000050 (2) — 000043 (2) +...). 
a) \, a 


15 


/ 3 6 
+ ft ors +-022715 (4) — 003120 (2) + -oor12i(4 
a \ a a ob 


~-000579 (2) “+4000355 ; 


~ 000241 (2). 
(a7 O 
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For the sectors whose angles are 90°, 180°, 360° the series are too slowly 
convergent for calculation, but the sums of the series are without much difficulty 
obtainable by the method used to obtain equation (46). 
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For y=90°: 
ee NO A Pine aay ih eH] 3 (4) 
ee (eer) =(" er) | 2+ (4+ (4) ) log (1+ (2) )= = 2 log,Y 4\a 


(2 
coszy| (7 1 vy") 1 
ty aes ¢ 4\a eee 


i 4) + (2) J2+1 Z 2 
GOES Ee 
+7 (24 (2) ‘ Pee neal Ata uly bn", (47) 


Instead of this complicated expression the series formula of equation (45) may be 
used directly for low values of ef Actually, in the calculations for the table given 
on p. 31, this method was used as a check on the accuracy of the work. 


For y=180", the case of a semi-cylinder flexed in the plane parallel to its plane 


side : 
ii Wan 7 1 fy\ Wa 5y¥ 2 eee aed LY +4) 
CH (npr) = 5-9 (4) + (7 BT) | ~22+5(1+% tan ee) (Z—tan Z) 
mpeg OL TCA aN Ra (48) 
For y= 270°: 


stb i) ti) Eel) Ct) 
thf (+2) PoC) 

“Yh 8(0) bern 
Ee) +N 


ven rfl) 8) IE 


and g=((X)°-(2 ie sosntcstcan onsen A) 
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For y=360°, the case of a complete cylinder with a radial cut perpendicular to 
its plane of flexure : 


Tf Wo Ney Wa\ 1 y Lig ee a y\~ 
Ke pr)>~? (at) eLJa\'~ala) ~ala) *(G) ) 


=P (é ay ) tan™ Z| 
a a a 


+(+ (2) 747 (2) = (4) ) tant /[¥1....- (60), 


In the cases of y=60°,y=90°, y=180", m takes the values 3, 2, 1 respectively 
and consequently there arise infinite terms in the series of equation (45). There 
are however corresponding trigonometrical infinities of the same order and the 


limits of the whole expression are, as they must be, finite. For y>180° the shear 


becomes infinite at t= . This is the case of the re-entrant angle referred to later 


on in this section. 


The accompanying Table III gives the collected numerical results and in Fig. 4 
are diagrams each representing the distribution of shear along the central radius 
for different sectorial angles. It would be impossible to show adequately in one 
diagram the comparative distributions of stress for the different sizes of the sectorial 
angle. It is to be remembered in considering the table and the diagrams that, in 
accordance with the axes of reference which we have taken, a shear of the material 
on the fixed side on material on the loaded side of the cross-section is positive when 
it acts opposite to the direction of load. 


Our diagrams provide the excess of actual shear, over mean shear in terms of 
mean shear. They indicate the point at which the shear is exactly equal to the mean 
shear and the value of the shear at the centroid. When y=0 we have, if A stand for 
area of section and W/A for mean shear, | 


sa TATA \ hee y\ 
fae |( W/A) 136 L412 (¥), 


or the limiting case is a parabola. The ratio runs from 0 to 2°56, or the excess on 
mean shear from —1 to 1°56. The curve is given in Fig. 4 (i) of this series. Almost 
at the half radius (‘508) the shear is equal to the mean shear, but on the right of this 
point, Le. towards the thick end, the excesses on the mean shear are considerably 
greater than the defects on the left, and we get light on why the thick end rotates 
upwards. We may look upon this curve as the parabola plus the vertical line through 
the point 1°0. This explains how the apex curve comes away as in Fig. 4 (ii) of the 
series, when y= 5", and the vertical line bends round so that we actually get negative 
central shears at the butt end in sector prisms of very small angle. A horizontal line 
through —1-0 of vertical scale indicates in all these diagrams the extent of the 
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Limiting case sectorial angle, Y=0° 
(Parabola of central shear) 
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Fig. 4. Diagrams showing the Distribution of Vertical Shear along the central radius of : 
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(xi) (xiii) 
a flexed beam having a complete circular sector cross-section, y being the sectorial angle. 
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positive and negative central radius shears. Passing to Fig. 4 (ii1) we find for y= 12° 
the defect on the mean shear still negative but very slight (—:023), the position of 
maximum shear is passing away from the butt end, and the position of mean shear is 
very slightly to the left of the mid-point of the radius. These features are still more 
marked in Fig. 4 (iv), y= 30°, where we have a considerable excess over the mean shear 
at the butt end, the position of mean shear is visibly separated from the mid-point 
and the shear at the butt end still increasing. The process is continued still further in 
Fig. 4 (v), y=45°, but the butt end shear in terms of mean shear has now reached a 
maximum and with Fig. 4 (vi), y= 60", has begun to fall. The shear at the centroid is 
now also falling; and the point of mean shear is moving rapidly towards the apex. 
With Fig. 4 (vii), y=72°, the tendency of the curve to flatten out horizontally is con- 
spicuous and this is accompanied by a fall in the butt end shear and the centroid shear _ 
with a still further motion of the point of mean shear towards the apex. In Fig. 4 (vii1) 
and (ix) for y=90° and y=120° respectively these features are emphasised. In 
Fig. 4 (x) for y=150° we are approaching a climax. The butt end shear is now less 
than the mean shear, the position of mean almost up to the apex and the flattening 
of the top of the curve now stretches across from 0°2 to 0°5 almost of the radius. 
A section between y=150° and y=180° would show the curve with a pointed apex ° 
like Fig. 4 (11) reversed right and left, and this enables us to interpret Fig. 4 (xi) which 
is almost a straight line. For y=180° the first point of mean shear has reached the 
apex, the zero shear has ceased and we have a finite excess over the mean shear at 
the apex ; the curve indeed should be considered as built-up of the nearly straight 
sloping line and the vertical through 0°0. The second point of mean shear has reached 
the neighbourhood of 0°8 and we have increasing defects on mean shear at the butt 
end. Fig. 4 (xii), y=270°, shows the revolution which follows on the re-entrant angle. 
The shear at the apex is now infinite; the shear at the butt end is now negative and 
probably approaches a maximum somewhere about 300°. The position of mean shear 
has also reached somewhat about its maximum in motion towards the right. In 
Fig. 4 (xiii), y= 360°, or the solid cylinder with a cut to the centre, the butt end 
shear.is again positive and less than the mean shear, the position of mean shear has 
returned towards the butt end, and the big shears are all towards the left, the exact 
opposite to what we noted in the sectors of small angle first discussed, but at the 
same time they are all to the right of the centroid, so that their total turning moment 
is big and the result a large couple turning the butt end upwards. 

The minimum torsional effect is shown in Table I to occur in a sector of about 
105°. This is roughly indicated in Fig. 4 (viii) and (ix) where the couple of the 
central radius shears at least, as measured by the moments of the areas above the 
horizontal line through — 1-0 about the centroid, appear mostly closely balanced. 

Of course the reader must bear in mind that these central radius shears are not 
the whole shearing system, but the central radius is the locus probably of maximum 
shear as in the case of the flexure of symmetrical systems. 
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§ 13. The torsion of the curtate sector for particular cases. 


The formulae for the curtate sector are very much more complicated than the 
formulae for the complete sector and require a large amount of additional arithmetic. 
This is however much lessened if the case of the complete sector for the angle in 
question has been worked out, for the calculation then consists of multiplying all the 
quantities concerned by appropriate factors. These factors, also, converge to unity 
with great rapidity if the ratio of the inner radius a, to the outer radius @ is not too 
large and if the value of y is fairly small so that m is large. This is the case, for 
instance, when y =12", which is given below. 


For the limiting case of a,—a, we can expand all the factors of the form a? — a,” 


= a 


in powers of e€= , and we obtain for limiting value 


0 


Wa Wa ; 

— See ? = ——— : : f = 1. 
7 = Fp (84°56 + 86°70) = Fae x 10624, if y= 
1°. A thin wedge-shaped curtate sector. 


The angle chosen was 12°, and the results of the calculations are given below in 


Wa . 
Table IV and graphically in Figs. 5 and-6. The column headed, ys, =7 / aT gives the 
(a—a,) . 


f W 
constants as directly calculated, and the following column headed ,=7 / “Tae ie 


deduced from the previous one. The first expression gives the ratio of total torsion (7/) 
‘ : ied ak oe s 
to deflection per unit length (5 = i“ ih dad it r 
height of section and #, the number in the second column of Table IV. It shows that 
for constant height of section and length of beam the torsion at first decreases slowly 


W,, where / is the maximum 
7 


Taste IV. The torsion per unit length in a sector of central angle 12° 
and external radius‘a for varying degrees of curtateness. 


Torsion per unit length 
Internal radius — a@ |———— oho =e 


Externalradius a Wa __ {W(a-4ay) 
ely 7) ver | EI 
| 0-0 2846 2846 
3 “2486 *3729 
‘6 | *2205 6614 
75) 2219 ‘8875 
8 2248 (toss 
9 2386 2°3861 
95 | 2686 53722 
‘98 “4528 23-6390 
1:00 106°24 aD 


i 


Norr. A positive torsion signifies a torsion in the opposite direction to the arrow in Fig. 2, p. 20. 
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and then, after a,=*9a, increases rapidly with the squatness. The second expression 


indicates that the same ratio may be represented by <5 x 3., where w, is the 


number in the third column of Table IV, and the values of ys, show that for constant 
breadth of section and length of beam the ratio of total torsion to deflection per unit 
length increases at first slowly and then very rapidly with the squatness. 


4 
Torsion measured on 
Breadth of Section 
3 
=> 
LIN 
sis 
S| 
Paps ok 8 9 
b 
1 
Fie, 5, 
O 
0'3 
SIN 
= ss 02 
aan 
i 
Torsion measured on 
011 Height of Section 
Fria. 6. 
00 


00 01 0:2 03 04 05 06 07 0:8 =F) .0'9 10 


Fic. 5. The torsion of a flexed beam which has as cross-section a thin (angle = 12°) wedge-shaped 
sector with varying curtateness, the torsion being measured on the external radius. 


Fie. 6. Torsion of the same beam as in Fig. 5, the torsion being measured on the 
central breadth, not the external radius. 
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Numerical example. 


As in the case of the non-curtate sector, we shall now consider the special example 
of a beam of steel 20 inches long having as cross-section a curtate sector of angle 12° 
and outer radius 1 inch (the outer are is therefore in all cases ‘209... inches). The 
values of the deflection and torsion when a load W is applied are given in Table V 


and the torsion per deflection for varying curtateness is shown in Fig. 6, the scale in 


degrees per inch being obtained from the scale of + lor there given by multiplying 


EI 
by the factor 430. 


Taste V. The deflection and torsion of a steel beam of 
curtate circular sector cross-section. 
Length of beam = 20 inches. 
Outer arc = 0-209 inches. 


Angle of sector = 12°. 
Stretch-modulus = 30,000,000 Ibs. per square inch. 


Internal radius _ ap be VAT tecaiot Torsion eda he 
External radius a =5 inches =rl degrees patie 
3 HI deflection 
| 
0 Wx °465 W x -0569° 0°°1223 
3 seed 1 «OUST. | 0°:1068 
6 x °580 x 0549° 0°:0948 
“75 x 681 x 0649° 0°-0954 
8 x “788 x 0762° 0°:0966 
o x, 1°353 <locome || 0°:1027 
95 x 2-509 x ‘2896° 0°-1154 
1:0 x 0 xX © 45°-68 


Norr. A positive torsion signifies a torsion in the opposite direction to the arrow in Fig. 2, p. 20. 

2°. A cylindrical tube with a longitudinal cut. 

If we take y=2z7 we get a sector which is a full circular ring except that it has 
a cut along one radius. The beam is then really a cylindrical tube with a cut along 
one side. This case, especially if the tube be thin (e.g. as below a,/a=-9), is of 
peculiar interest as it represents the extreme of asymmetry in a circular sector. The 
method of approximation for a,>a mentioned at thé beginning of this section having 
failed, we have to use direct calculation and, both y and the ratio “t being large, the 
calculation of the series is extremely heavy. The factor of 7 in equation (39) turns 
out to be the difference of two nearly equal numbers and to get three figure accuracy 
in the final result it was found necessary to calculate a large number of terms of the 
series to 10 decimal places. 

If we take a,="9a we get 

Wa Wa 


= —— ? S Se : hg ele 
T= BL * 2206 (1+ 1) ay 850°75, for 4 =4, 
_ Wa 


as compared with | Rockey oy x 1°057 
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for a non-curtate sector of the same angle. This is an extraordinarily large value for 
7, but it is well known that a hollow rod loses much of its power to resist torsion 
when it has a cut down one side and it might therefore be anticipated that the couple 
due to the great asymmetry of the flexural shears would bring about a correspondingly 
large torsion. 

We can illustrate this result by taking the value of 7 for e=(a—a,)/a, very small. 
The deduction of this value is of much interest as it involves the summation of a number 
of noteworthy series which, however, cannot be considered here ; ultimately we find : 


Wa 4:7124(1+) Wa _ Wa 5°8905 


ne 9 é ER EE BC ETT 
indicating the high value of 7 as the thickness diminishes indefinitely*. 


Experimental Verification. 


With such a large value of 7 it is easily, possible to see the twisting of such 
a beam under flexure, and some rough preliminary experiments were made with a 
piece of piping used as a casing to electric wiring. This is manufactured with 
the required slit. One end was fixed in the chuck of a turning lathe, a solid 
boss having been put inside the tube at that end to prevent buckling. The load was 
applied by a stirrup arrangement fastened to the other end. The deflection was 
measured by taking the vertical heights of the bottom of tube from the lathe-bed or 
from the slide-rest. The torsion of the load-end relative to the fixed end was measured 
with the aid of a light pointer fastened rigidly at the load-end moving over a plate 
fastened rigidly to the fixed end. The effective radius of the pointer was 10 inches 
and the torsion was measured on a 12” length of the 15” pipe so that with the largest 
observed torsion of '159 radians or 9° the pointer moved through 1°59 inches, and the 
real torsion was 3. 159=:'199 radians. 

The dimensions of the tube were: Outer radius = 0°50 inches, inner radius = 0°45 
inches, length (excluding the part used in fixing) = 15 inches. 


TABLE VI. KHapervments with a piece of split piping. 
Deflection and torsion with various loads. 


Experimental Theoretical (H =27,000,000) 
: : 
Sha w Average deflec- | Average torsion=7l Deflection in Torsion in radians 
tion at end in -- ——] inches (linear - 
inches=6 In radians In degrees |andecubicterms)| q,—+45 | a)=-425 
0 000 000 Oe 000 ‘000 “000 
37 048 O15 * 52! 037 041 O17 
13°7 ‘131 059 py sy! 138 ‘151 063 
23°7 215 105 Gaels "239 262 108 
33°7 “305 152 8° 43’ 339 372 153 
43:7 387 show 11° 24’ “440 488 Og 


Wi 
* Tf ¢« be put 0:1 we,have 7 = G7 589°05, which is of the same order as the true value above, but 


indicates that ;/5 is not small enough for us to neglect the next term in «. 
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In the first experiment the deflection and torsion were measured only at the end 
of the beam. The results are given in Table VI and represented in Fig. 7. 


Torsion in radians, 7/ 


Deflection in inches, 6 


0 5 10 15 20 25 380 35 40 45 
Load in lbs., W 
Fie, 7. Deflection and Torsion with varying load of a split pipe*. 


Deflection in inches, 6 


Distance along beam in inches, z 


Fic, 8. Curve of Flexure for split pipe. 


It is evident that both the deflection and the torsion are proportional to W, but 
there is a great divergence between the observed values of 8 and the values that 


are given by direct use of the Eulerian expression UG for 6, # being taken 
to be 27,000,000}; Ibs. per square inch. This leads us to question whether the linear 
term s,z in the expression for the displacement is negligible as it is often taken to 
be in beams of moderate length. As we shall see the theory does not give much 


* The deflections were taken on a 15”, the torsions on a 12” length, so that the latter need to be 
multiplied by 1°25. 


7 This choice is based on calculations given later on in this section. 
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help because the nature of the clamping is so uncertain, but an experimental test is 
easy and obvious. If we measure the deflection at various points along the beam we 
ought to be able to pick out the linear term if it be not negligibly small. Observations 
were therefore taken of the deflection of the beam with a load of 43:7 lbs. at three 
points along the beam. The means of several trials are as follows. 


TaBLe VII. Eaperiments with a piece of split piping. 
Observations to determine curve of flexure. 


Distance along beam (z) in inches: 0 44 9 16 
Average deflection 6 in inches ...... 0 ‘070 166 *313* 
Approximation of equation (51)...... 0 068 169 313 


These are plotted in Fig. 8 on p. 39. If we assume the theoretical formula 


8, = 82+ ue eae 2) a aa iS 3 (51), 
and substitute the values of W, J, I, we find by the method of least squares that 
s, =°014, 
and E=27,000,000, 
giving 6, = 0142 + 000722 — 0000162". 


This value of # is quite near what might be expected for this grade of mild steel. 
The large value of s, appears at first extraordinary but it is even less than some 
theoretical forms of clamping provide; with it we should have to take a length 
of over 65 inches before the linear term was less than 10 °/, of the cubic term in 
the expression for the end-displacement. Exactly the same experiment was made 
on the pipe rotated axially through 90°, so that 
instead of having the cut to the side as in (a) of 
Fig. 9 it was vertically upwards as in (b), and 
the least square approximation determined as 
before. This gave 
= "0082, 
E = 27,000,000, 

the same value of £ as obtained before. Rougher (ay. (5) 
observations for other positions of the cut showed, Fia. 9. 
as one might expect, that the s, is a maximum 
when the cut is at the side and a minimum when the cut is at the top or bottom, 
these positions giving roughly the same deflection. | 

* This small value as compared with that in the previous experiment is, no doubt, due to elastic 
fatigue. In the first experiment the load was applied by additions of 10 Ibs. and between each a 


considerable time elapsed while the measurements were being made. In the second experiment the 
43-7 lbs. load was put on and off at short intervals. 
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Using the value thus obtained for H we get as our theoretical value for 71, 
when W = 43°7 lbs., | Tea, 


Thus the theoretical values of the torsion are more than double the observed. This 
may perhaps be accounted for by the roughness of the observations, especially of 
the thickness of the pipe, for it is evident that, in the neighbourhood of a,=‘9a, 


Wa .. : yok : 
the value of 7 ET will change with extreme rapidity. For example, if a,=°425, 


Wa : : 
t/ 37 = 214°16 + 214°18y = 267°70 for »=4 and the torsions are then those in the 


last column of Table VI on p. 38 agreeing almost exactly with the observed. Moreover 
a tube of this form made out of rolled plate is probably fairly aeolotropic, and this, 
_as will be found, changes the value of 7 very largely. There is further the difficulty 
of arranging the load so that it acts through the centre of gravity, and above all the 
essential difference of the methods of fixing in practice and in theory. 


§14. The nature of the fixing at the end of the beam. 

The experimental results of the preceding paragraph lead us to consider more 
particularly the nature of the end-conditions—what parts can be fixed and how, so as 
most nearly to resemble the actual conditions of practice. Physically we see that any 
extended fixing is impossible with our solution, for, as shown in § 9, it is at the fixed 
end that the lateral displacements are greatest and to them must be superadded 
the longitudinal distortion which is the same for every cross-section of the beam. 
Mathematically we see the impossibility of extended fixing when we consider the 
fewness of the arbitrary constants which arise in our solution and which we alone 
have at our disposal to describe the manner of fixing. Bearing our limitations in 
mind, we proceed to treat some particular cases of theoretical fixings. 

From the expressions for the displacements [§ 4, Equation (20)] we derive for the 
end cross-section, where z= 0, 


Ww 
U=sr aml (w—y)—yy ta! 

W j \ | ; 
v= az neyl+ yet B Sea Alec coneoea sty: (52), 
w=yx (a, y) — py oy — Blat aly ty! 


and, arranging the differential coefficients in a Jacobian scheme, 


1 CeCe an W W ER aa Dar 
dn’ dy? dz = pr — Ry y = TY — F790 (H eat uate 
dv ov dv W ] Wey De ee 

Bx? dy? dz EIN TY Bie eS Pas 

ow ow ow Xe Py OX pe ye srk be ae 
e/a fe Sa a 2 i ae yd hoes 


By means of the equations in (52) and (53) we can fix certain points and certain 
differential coefticients at these or at other points. 
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(i) The fixing at one point only, in particular the centroid of the end cross-section. 

Suppose we aim at fixing as far as possible the neighbourhood of the centroid* 

of the end cross-section, namely (0, y, 0). Firstly, to fix the centroid itself, we 

have from (52), Wa WwW, PaaS 
= pp oyl—-yyte= 


v =f’ =0 A SANE AOR atid aa (54), 
w=x(0, 9) +aly +y' =0 
Ou oO 0 : ; . ns 
and automatically = ay = = =0 for any point on the plane «=0, which is the 


neutral plane. 


; ‘ ov ow ; 
To determine a we may assume either a, OF ay) to be zero, since we know that 
a=a" for our case of symmetry about the axis of y. This makes a=a’=0, and 
means that there is to be no rotational displacement of the beam as a whole about 


the axis of a. 


: ou ov , Wien : 
To determine y we may take By =e, making y= — R ply: This really 
determines that there is to be no rotational displacement of the beam as a whole 


about its longitudinal centroidal line. 
From equations (54) it follows immediately that a’= B’=y'=0. 
The only remaining constants are 8’ and B”, between which there is a relation 


given by either of equations (37) and (38), so that fixing the one automatically fixes 


(3 
the other. We may either make ze =(0a r= 0. 


Ls ou 0, Apart from the mechanical difficulty of introducing a means of 


fixing which will fix an element in the direction of z, this would make the linear 
term (sz) in the equation of the curve of flexure vanish, and this is contrary to our 
experiments. We shall therefore make use of the second alternative which was that 


adopted by Saint- Venant. 
y Ve 
2. Lees or ) — A aby 
y 


Ox 0x LI 
This is the same as ay Ww Bus, 
v00), EI? . 
: 0 ; ; et 
since —=_= 5 on the radius of symmetry, or, substituting our value of x from 


equation (30), 
os B 

| = (ar ee a m cos m@ —(B — 8B”) cos de ree 2 cos 20 
i ee f 


oe 
Ae E cos $y * 
any | $+(4—-4n) —4 aE ; er ?— B” =0, 
3hT| 4" (4-3) 1 08 at aad B 


ay =) 
cos $y 


* No property peculiar to the centroid being used, the equations which follow hold good for any 
point whose coordinates are (0, 7, 0) on the central axis of symmetry. 
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or 
mA B, TY W cos ty] _ W _ 
A qpme-2 = eee ate jin oY Ape Wie oe 
Ex nf eT i ree 1 a|? +(4—2) COs | Aa Ris ME 
giving 
bee Tab Be TY W cos $y] -, 
B a e (4.9 on a m+ cos y 7 En 4 t+ a zn) ae Yo sesisones (55). 


Collecting our assumptions and results, we have : 
If the centroid of the end cross-section be fixed so that 


= Via a= () 


on 2 _ wn _ 
ie Oa Pe ea Sa ey ea oe (56), 
m0 _2w_ u_%_. aw_, 
ery ame Outs. ony kV Ont 


then the displacement at any point (x, y, z) in the beam is given by 


W ae ; Wire 
w= ryt ayy | 49 (ey) (1-2) +e —4 | + rn uy 


= 7jmn1 ep Ty a 1 
+2 B (4, y ea "+ Cosy = eo - ( ene 


Ss! 


cos 
(£=47) - 1) 4 | 
COs 
IY SN 57), 


V = Tet + ney (1-2) — py nly 


V I,2 2 Ox ie 


If we consider in particular the curve formed by the line of centroids, that is the 
curve obtained by putting «=0, y=y in the equations of (57), we obtain after a 
little reduction 


w= x (x,y) — 


where tab is the vertical shearing stress at the centroid of any cross- 
section (from equation 45), and 


ard WwW Wa?zl 2 B yi 
U(Y) = yy (2 —$2)+ hides is (409 re ok =| m | 
Wa ¥ y ty cos 37)|2 | 
ra «lay ees anaes | 
assuming uni-constant isotropy, (58 

W i vie | Rb EDS 
=A ($l —Ad) + is Z | 
| 


v (y) =0, w (y)=0. 
(a) Case of the complete sector. 
For this case the distribution® of the vertical shear along the radius of symmetry 
has been worked out, and we can determine the values of za at the centroid or at any 


* See § 12, Table IIT and Fig. 4. 
6—2 
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other point with sufficient accuracy by interpolation. We thus have the following 
table for the value of the coefficient of the linear term in 1 (4). 


TasiLeE VIII. 
Angle of sector: 5° 12° 30° 45° 60° 72° 90° 120° 180° 270° 360° 
saad ie ube x |°00115 | 00674 | 04144 | 08896 | +15127 | -20726 | -30371 | -49983 | 89385 | 1:20921| 
Bb | : . 
| | 


It might be anticipated that a better approximation to what would arise in 
practical fixings, when the sector has a small angle, would be obtained by taking 
the theoretical point of fixing to lie between the centroid and the circumference 


: : 20 
end of the radius of symmetry. The nature of the corresponding values of — can be 


pe 
immediately seen from Table III. Such a point of fixing would enable us to avoid 
the infinite linear term in the deflection for the split cylinder. 


(b) Case of the sector, complete or curtate, when y = 2r. 


In this case 7 is at the origin and when ¥=0, 2a; becomes infinite for the com- 
plete sector. Moreover we must always have the relation of (38) satisfied whatever 
be our theoretical method of fixing and it is clear that for our case of y= 27 the value 
of B— B” will have to be infinite, since 7 is not zero. 


These mathematical infinities need not deter us from applying the theory to 
physical problems. In the actual experiments we cannot get angles of exactly 27 and 
we have no means of practical fixing in accordance with the theoretical methods. In 
fact the true interpretation of these infinities is that we need not be surprised at 
relatively large numbers arising at the corresponding points in physical problems and 
it was the large value found in the experiments with the cut pipe which led us 
to the investigations in this section. 


(ii) Laing at two points—case of cylindrical split pipe. 

Instead of fixing one point and certain directions through that point we may use 
the arbitrary constants in the solution to fix two points of the terminal face. The 
choice of the two points would naturally vary with different kinds of sector and we 
shall devote our attention immediately to the cylindrical split pipe. 


We shall choose the two points (+4, 0, 0), ie. the two points where the plane of 
flexure crosses the outer circle of the end-section. .The equations for u=v=w=0 
will then be 

ed AG Ps 
U= panel ge 


v= +ya oa 0 Ee (59). 


r=a ee 
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From these we deduce immediately 


/ 


a= LT ancl [TUG eects ey Sed 
and B’a=x (a, 5) , 


: OL ; ‘ 
To determine a we may put y= ° at (+a, 0, 0)—making a=0 and virtually 
preventing a rotation about the diameter through (+a, 0, 0). 


The equation for B” is from equation (30) 
= Bin \e 1 OUT, ae ' 
Ba = mh (4 mo at: =| sl —(B- B") a — 3KI [# a3 (F a $n) | a, 


= Pie W 2 
whence B= oe (4,0 a 4.) sin (27+ ae 1 — 37 (i —B9) Od. von ewes (60), 


or, substituting for A,, and B,, from (34), 


1 9 (a) (=) 
& Wet ly 8 x 520°6 (1 +7) a]. a 


“| a* KI Qar« na MW (m? — 4) i. +a sy" 
O 
m+ 2m 
B42 3—6 1-2 (4) le, fs 
a ae ey ein CE Ve 
m (m? = ‘jam m n (mi —9) Ns as in 
a) 

1W : 

—3 ay —tyn)a@ Coe eer are sorcerer eserneeserneeeesees (61). 


With these determinations of the arbitrary constants the displacement equations 


become 


lag W 
u= — Tye + elon (2 —y') (1-2) +312 —42']+ Be— a7 bya 


W : 
v= Tea + Fre acy (1 —2) ...(62), 


W HN 7 
w=x (ty) — pyle (eh?) +2y']- = x (a, 3) 


and in particular the curve of centroids, which for the sector in which y= 27 is the 


same as the curve of centres, is 


CS eats 
Uy = arr le? — G2 2) + peal OT AE ey ore (63). 


The last constant term is of course the result of the anticlastic curvature and as we 
shall show is of no importance compared with the value of B/. The value of £ is 
identical with s, as in all cases where the centroid coincides with the centre. 


Returning, therefore, to our value of 8 in equation (61) we find that the first 
ten terms of the series are needful to give adequate accuracy when a,= 454, but the 
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only terms that are of real importance are those containing the large constant 
520°6(1+ 7). The result of the calculation is that 


B=s= "po x 907-05 Dee aiee acheter (64), 


an extraordinary value when compared with the values (for angles other than 360°) 
in Table VIII, and with the comparable values for the case of the beam of rectangular 
cross-section calculated by Saint-Venant™*. 


This large value is; however, confirmed by our experiment described in § 13. For 
when we calculate 6 for a pipe of the dimensions used in the experiment, we find 
that B= "0240. This is somewhat larger than the observed value which was ‘014. 
However it is of the same order when we compare it with other cases. 


It is interesting, for instance, to compare it with the value of 6 deduced by the 
same method of fixing if the pipe be closed. It is easy to show for this case of the 
complete hollow circular cylindrical pipe that 


UA Gestation. ao, vee Wawa ye 
x= gy +30) ((@ hag) Se) sin O+ spar sin 30—(8—6”)rsin@ ...(65), 


the notation being the same as used for the curtate cylindrical sector, and that the 


value of B, when the points r=a, 0= ee z=0 are fixed, is given b 
P 5) g yy, 
Wis on pipes 
B=77 (84 ot Og Moamop caw + Muduewieep ase: teeaee (66), 


or, when @=0°50 inches and a@,=0°45 inches, 
B=-000053, 


about 345th of the corresponding 8 when there is a cut in the pipe. 


Now the plugging of a split pipe and the clamping of the plugged end in the 
chuck of a lathe may not very closely correspond to fixing indefinitely small elements at 
the top and bottom of the external circumference at the fixed end of the pipe. It 
is worth while inquiring what value 8 would take, if we fixed external elements on 
the horizontal diameter at the clamped end. We may take as our conditions: 


u, v, w=0 for z=x2=0, r=a, 0=0, 


and again for r=a,0=7. 

WwW 
From (35) we have U=0= FF (—4dyv) ltyata’, 
whence flow y=0, a =4 sis nl, v=0, B’=0; 


w=0= x (a, 0)+a"at+y, 
w=0=y (a, 7)—-a"at+y’. 


* “Memoir on Flexure,” § 29 (Journal de Liowville, Series 11, Tome 1 (1856)). Quoted in Todhunter 
and Pearson’s History, Vol. 11, Section [96]. 
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Now x (a, 0) is zero, and 


B 
x (a, 7) = 5 Nee a” + A) sin (2i+1)= 5 
does not involve any additional unknown constant. 

Hence these equations determine 


" 
a 


1 
= 3a X (a, 7) and y/ = —4y (a, 7). 


It still remains to select a condition which provides 8. A suitable limitation seems 
to be that dw/dax should be zero at non-split edge, or 
we (94), 2 gh 
dx \rd0/gu0 
: s Wa 
which leads to Bas 2 (4,0" te a) m+ra—sn se : 


20°6(1+ 7) for 7, we have 


Substituting for d4,, and B,, from (34) and the v ie 7 


; yy (ay (<2)" 
jae ( Qa we A i -(“)" 


ah m+3 a} \ 2m 
sa Di 22 eo ies 
ah ah 


ops! [Q, ay 


\a@ 


+ 520°6 (1+7)—43n}» 


il 
m—1 Werte 


or, if »=4, p=" n" x 959 675, 


a result comparing fairly with the previous high value 927:931 obtained by fixing 
elements on the vertical diameter. It seems necessary to infer that with perfect 
clamping we should reach a value of 8 comparable with that given in (64). 


Experimental verifications of the foregoing theory. 

It will be of advantage at this stage to summarise some of the results of experi- 
ments made with a view to increasing the faith in the theory. In all the cases 
considered here, the deflection (6,) was observed at various points (z) along the lowest 
generator™ of the cylinder, and the constants s, and s, in the formula 


Ol Ope Sie ae ee Vet ches eu dn cves<ceuev-W bay ses (67) 


were deduced from the measurements by the method of least squares. The results 
are given in Table IX. Remembering that s, seems naturally to arise in the form 


8, = (Wa’/ET) x a quantity depending on the shape of the cross-section, 


2 
* s, will therefore differ from 8 by a term of the form a ty. This term would have no sensible 


effect in the case of the split pipes, but’ would mean that s, is less than 8 by about 6°/, in the case of the 
non-split tube. 
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where @ is half the typical vertical dimension of the cross-section of the beam*, we 
shall write, after Saint-Venant, the total deflection as 


Whe La Ba 
o= sar! (1+ Fro.) midle s:6)<.xiselbiolelie 606.9 v's epee sap (68). 


The quantity o, is tabulated as well as the ratio s,/s,. For purposes of com- 


parison, of course, a, 1s superior to s,/s;. 


TABLE IX. Haperimental results on the nature of the deflection 
with various beams: 6,=s,2+ 8, (442-42). 


Description and dimensions of beam - 81 _ y=" ate CE 
Split pipe (soft steel): a =0-°50 in. 
cut at side | %=0°45in. | 0136 000096 140 560 27,000,000 
Z=15in. (0240) 
Same pipe: cut on top a =0°50in. | 
| a =0'45 in. | 0082 000096 | 85 | 340 | 27,000,000” 
Z=15 in. (0068) 
Soft steel pipe without a =2in. | 
cut @) iy = zy in. 00181 000236 | 7-7} 54 |- 23,000,000 
1 =9:5in. (-000061) 
Bass-wood cylinder gt se 0151 | -0000681 | 220 | 570 | 1,470,000 
with saw-cut to t =382 in. 
; ; (00046) | 
central axis | 


The experiments were only rough and little was expected of them beyond some 
slight verification of the theory. The values, however, of the Young’s Modulus, £, 
are in accordance with what we should anticipate for the materials used. The chief 
divergence lies in the values of s, when compared with the theoretical values of 6 in 
the case of the solid pipe and the wooden cylinder with the saw-cut. But the nature 
of clamping, by a chuck for the uncut pipe and a simple clamp for the bass-wood | 
cylinder, is really not that of two fixed points and may well provide somewhat more 
closely the high values of a central fixing: see (b), p. 44. 


We may note: 1°. the difference between the o,’s for the split pipe when the 
cut is at the side and when it is on the top (the arrows in the drawings of the cross- 
sections being in the direction of the flexure) ; 


2°. the small value of o, for the complete tube although it is not so small as the 
method of fixing by the two points (+4, 0, 0) described above would provide ; 


3°, the o, for the wooden beam with the saw-cut up to the central axis. This 
is the physical approximation to the non-curtate sector of angle 27. The o, is large 
as we should expect from a consideration of the order of the ratio of stretch-modulus 
(£) to slide-modulus (4) for wood. This is considered more fully in the next section of 
the paper. 


* In our case the outer radius of the sector, in the case of a beam of rectangular section the semi- 
thickness, and for an elliptic section the semi-axis parallel to the direction of flexure. 
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IV. SoLuTioN FOR CERTAIN AKOLOTROPIC CASES. 


§ 15. <A solution for aeolotropic material. 


Reverting to our original equations in § 3 for a uniformly aeolotropic substance, 
we see that the main change in the analysis that a consideration of aeolotropy 
requires will be introduced by the equation (15) of p. 10 


instead of the usual Laplace equation. It is natural to try a projection, which will 
transform this equation into Laplace’s equation. Let us project, then, with the 


substitution 
om, [He EE he it Ee REPS Ay POE (69 b). 


This projection depends on the nature of the aeolotropy, circles being changed 
into ellipses of semi-axes proportional to /, and ./#,, but in such a way that there 
is no change in the lateral dimension of the beam. 


If the boundary be f (x, y) = 0, 
the transformed boundary will be 


b(t rat y')=9, 
or ' (x, y') =0, 


if we agree to denote everything in connection with the transformed problem with 


dashed letters. Also 
Ad J/* op 
cos (a’v’) = Tey) a Tey IG: | 
Ue SEG) +] 
p 


a 
=, |" cos (x) x x — eee PON OE Oo Ae (70), 


py \O a2 
and, similarly, / ( a) “it ¢ fy 


Ox 


s (y/v") = Cos (yr) X — | vee cette cee ee cece neces Walp 
cos (y cos (y Je ay + (3) 


Thus the problem of the aeolotropic case, namely : 


To find y (x, y) and 7 to satisfy 
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with the boundary condition 
by cos (av) + He cos (yv) =7 (wy Cos (av) — pw cos (yv)) 
W =. ER Se ke een eye 
aa (Same a 1 a abla) He r) cos (av) + (H—7,p) xy cos (y) | 
— p, (B—B") cos (av) + py (a — a") COS (YV)  seceecreeeecneeeeceeeeeeeeees (72), 
and the Conditions I and II, 


=|x C08 (YP) dS <.0s\ocateransnkenehemmee eee (73), 


Baas 2 2 ox W 3 2 
-w9=|||- (x + mg) + (aoe oy - myx) By E ~ Bm) BY 


Wi En 
MF ia 2 pe 


is projected into the problem : 


2772s) Ba YU — pr (a— a’’) HL — phy (B—B") y| dady ...(74), 


To find y’ (w’, y’) and 7’ to satisfy 
Oo” O° 
ON anaes 
da’? * aa yl!” 
with the boundary condition 


/ 0 / 
as cos (ov!) +55 cos (y/v’) = As (y' cos (a’v’) — ax’ cos (y'v’)) 


a Bi Hs | (Bee ry, B= Mba — 2M bs *) on , B= hE 
4 ef COB LOY) ae 
EI \ 2 2 py q aes 


ip ar, cos (a7 v') +a = a)! GOSay u)) Vi hin veeseacke peste eeeeeeen (7a); 
and the Conditions I and II, 


0 =x CORAM BIS. tees. batssas seer omen e476), 
Aes, , OX ! 
a wy=|{[ om (ol + yf) Spa (a teh 2) 
W 3 Ww (E—y Waal 27). M2) Pats 
— py (E—¥mm) 7 aly! a oh 2 ps a 
acta oun ot eae een y'| /e ded dif ssc ae (77). 


This projected problem is analytically exactly the same as for the isotropic case™. 
The only difference from the latter lies in the different constants in the three condition- 
equations. Hence any solution that has been obtained for an isotropic case can be 
used to obtain a solution for a connected problem in an aeolotropic material. There 
is however the limitation which must be emphasised (to the depreciation of the 
method), that the form of the soluble cases depends essentially on the solubility 


* Note that Z=(p,/j1.)? /’ has been retained and is therefore = + (4/42)? (a* — a;*) (y' —sin y’) for the 
‘curtate elliptical sector.’ 
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of a corresponding isotropic case and also on the nature of the aeolotropy—to be 
precise, on the ratio of the two slide-moduli p,, »,. The result is that the soluble 
cases cannot always be for cross-sections of common occurrence, but must be for 
distortions of such. In the uncertainty of the elastic constants of aeolotropy, how- 
ever, this disadvantage may in many cases be of little importance. 


$16. Case of the curtate sector cross-section. 
Let us take the projected figure to be the curtate circular sector which is 
enclosed between the circles 
gl? 4 1 yl” = (2 
eA afm Oe Wish getacen a erences domes (78). 
and the lines y = +a’ tan ty’ 
Then on projecting back by the transformation 


the cross-section given will be a ‘ curtate elliptical sector’ enclosed between the ellipses 


P2 0 2 2 
—x+y=a 
Py 
CS oh 4. ae on Rate alas (79). 
Py 
and the lines oe +, [8 « tan ty = +a tan fy, say 
Fa 


Following the method of treatment used for the isotropic case in § 6, we get the 
following solution for y’ in the projected problem : 


Ante a a abe , 
1 alm m / / 1 y 1 on DAN 
x= a 5 (4 m! pay 72) sin mo’ — —(B- B i [2 7’ sin 16 +7! fe anes , sin 


ay peo 

+ : a ak eS ~ BPs =F Mth +8 pi 4 EW El a | 7? gin 30’..... .(80), 
be Be COS BY 

where 


7 m+2 m+2 

ee (otk art a 

mie OF 2 / 2m 2m 
Pf. (m’>—4) my! a —a, 


| 
W si cos 4,/ 3 h— No pes 1 3h — 97. My il qm gh aa | 

“t, ( i 1) oF . nae BAe ae ate Tage! | 
| 

| 

\. 

| 

| 

| 


ey. Be my’ ite m—1 ie a G 
B’ <; 87! py ( a 1)’ eg re Bir a"~*) 
c. Pe (m* =a 4) a Vg fete a a” 


al (e2i): cos ay 3H — 1 fs iS _ 8B -99. ps 1 ) 
i my’ fp m—1 i m —9 


. +3 +3 mM—-3 m-3 
a ag" *? (ar aos ) 


2m 
ea 


and m= (20+ 1) %, goa Te inte 
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Condition I will as before (p. 17) be identically satisfied for symmetry perpen- 
dicular to the plane of loading applying Condition II, and remembering the value 
of I (see ftn. p. 50), we get as equation for the torsion 7’, 


0 64 (ae? as ays 
/ 4_og4 I_eoto / - : 
T (a a7 ) (y tan of ) ae 2 (m ny 4) (m ae 2) my! a” ei oR 


4+=0 


¥ as 64 a! ig (ans — ny] 
i-o (m* — 4) (m—2) my’ a” —a’™ 

; Sony (Eh 

=o (m+ 2) my’ \ fs m—1 fs m—9 


es Wr 


e (am** ree Gate) een =a 


2m 2m 
_— a 


a 


_§ _2cos gy (2 =A, 1 _ 6H~ 36M, 1 
=o (m— 2) my’ [hs mL fs m—9 


ye be! (a™*—a,”"*) (a™* —a,""*) 


2 H—6y.p, . 
Beene si LHe by’ tan 3y! (a —a,’) 


2 E € sin x) pee Ly/ (a? — t,°) (as | 


 — Ob, 
Now projecting back, so far as this equation is concerned, consists merely in 

‘ malig "sad 2 Wee : ! Vie =4 Ma aI retaini 

replacing* 7’ by 7, and y’ by y, where y’=2 tan a fe tan by) and retaining J, so 


that with these changes the equation (82) gives the value of the torsion accom- 


panying the flexure of a beam having the curtate elliptical sector cross-section given 
by the equations (79). 


§ 17. Solution in the case of transverse isotropy. 


Transverse isotropy consists in the slide-modulus and the Poisson ratio being 
independent of the direction perpendicular to the longitudinal axis of the beam. We 
may write p,=p,=p and 7,=7,=7, but we do not necessarily have H=2p (1+7) 
which would give the case of complete isotropy. 


Applying the method of the previous section we see that the projected figure is 
the same as the original figure in this case, but we can deduce the result from the 
aeolotropic case and have immediately from equation (82), 


* Of course ‘replacing’ 7’ by +, and retaining J does not involve equality of + and 7 or of J’ 


and J. 
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’ fo 0) 64 (a™*2—@,™**) 
T| (a'—a, —ta ree? te 0 
( ) (y an Y) a oe (m* 7 4) (m ne 2) my on” — Cb, 2m 
rp S tem 64 , a hy ‘(a m2 Gea): 
iso (m? — 4) (m—2) my a" — a," 


Pet Ona? Cogs bE 1 ee 1 
Ast Ss sy ui wy, 
re 7 cd 6 = ee yp —— = 
fl e (i + 2) my (( be ~ 4) m* — 1 (6 BK 36n) mn? — ‘| 


(Ons = ee) cern Oy ral 
x “YO : Z2uT8 


-y 2 2m 
a" — a, 


2 cos 4 iy 
Se 27 ((6 E ee tn) artes pail as BY 36y RAS 
i=0 (mm —2) my mn m —1 nm m—9 


e a Oy ate dat aes *) (oie me Gases 


E Ce =e ar 
l - 
15 15 te in -3n) cos $y tan ay (a —a,°) 
2/1E 
er 5 ime 3m) sin gy (4° — a’) 
25 sin oe (a* —a,2) (a — a8) 
Tae =y aih eae fet sO : ee 
3 (1 y )sin 27 ees | Se ee ee i rr (83). 


A comparison with (82) shows us that this result is absolutely the same if the 
appropriate H/y, and y, be used. In other words the difference of the two cases lies 
solely in the distortion of the auxiliary cross-section. 


§ 18. Note on the torsion of an acolotropic beam with flecure. Cf. § 9 
The equation for the vertical displacement at the point (a, y, z) in the aeolotropic 
case is (equation 14) 


W 
U= — Tye + ard (1-2) (ma — gy’) + le —§2]—yy + Beta’, 

the only change from the isotropic case being the introduction of the two Poisson’s 

ratios y, and y,. When we consider the displacements on the horizontal axis of 


symmetry the only Poisson’s ratio occurring is the y,, and we have as the torsion in 
the case of the beam discussed in the previous section - 


— b= —Tl+ We Ce Os) et Jens ae aeaelin ea Meds (84). 
Compare equation (42). 
Thus the torsion due to the anticlastic curvature depends only on that Poisson’s 
ratio which is the contraction in the direction of y (i.e. perpendicular to the plane of 
flexure) due to an extension in the direction of z, the longitudinal axis. 


§ 19. Applications to certain aeolotropic materials and experimental results. 

The interest in the aeolotropic case does not lie so much in the immediate 
usefulness of the formulae as in the estimation of the order of magnitude of the 
changes that the aeolotropy of a commonly occurring kind is likely to produce. 
Aeolotropic materials other than those of regular crystalline structure are by nature 
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so variable in their aeolotropy that it is impossible to give dependable values to their 
elastic constants. In the case of the common woods for instance all depends on the 
part of the tree trunk that is selected, and in any one plank the aeolotropy may be so 
variable as to upset even the nature of the mathematical result proved for uniform 
aeolotropy. Thus in an experiment that was made with a plank about four inches 
broad and bevelled so as to have for cross-section approximately a circular sector of 
small angle, it was found that the torsion caused by flexure was most pronounced but 
in the direction contrary to that given by theory. The annual rings in the wood were 
noticed to be much wider apart towards one edge of the plank and, this explanation 
having suggested itself, the experiment was repeated with the ‘complementary’ piece 
of wood, i.e. with a plank in which the rings broadened towards the other edge of the 
bevelling. This time the torsion was of the same sign as in the mathematical theory 
but the magnitude was quite out of proportion to what might be expected with 
uniform aeolotropy. 


It is conceivable that this property of the variability of the elastic constants in 
wood might be made use of in those cases of practice where it is important that the 
flexure of a beam of asymmetrical cross-section should cause little or no torsion or a 
torsion in a certain direction. Although it might be difficult to select single planks 
that would give the suitable torsion, it might be easier to make compound torsion- 
free beams of several layers of different pieces of wood. 


Asa rough guide to the magnitude of the change due to aeolotropy, it is easy 


to consider the case of transverse isotropy (§ 17). If f is fairly large so that we can 
v} : 
neglect multiples of 7 in comparison with —, we see from equation (83) that the 


value of 7 is roughly proportional to ae The value of : for an isotropic material is 33 
but for wood it is much higher, if we take H in the direction of the wood grain. In 
the case of beechwood, for example, we can estimate* that ta is about 6, and for 
pinewood about 30. ‘ 


Some measurements were made of the torsion produced in a round beam of bass- | 
wood having a saw-cut through to the central axis. The beam had then approxi- 
mately the cross-section of a complete circular sector of angle 360°. The radius of 
the beam was 2 of an inch and the length of the flexed part about 36 inches. The 
flexing weight which was used varied up to 12 lbs. 


To measure the amount of torsion one end of a stiff wire was inserted in the saw- 
cut at the free end of the beam so that the wire stuck out horizontally and at right 


* The estimation is based on the values of the stretch-modulus given by Wertheim and Chevandier 
(Comptes Rendus, T. xxi, p. 663, 1846) and of the Poisson’s ratios found by Mallock (Proc. Roy. Soc., 
Vol, xxix, p. 157, 1879), and is made on the assumption of the ‘partial ellipsoidal’ elasticity of 
Saint-Venant (cf. Todhunter and Pearson’s History, Vol. II, Part 1, § [314)). 
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angles to the beam’s axis when the beam was in its normal position. The rotation of 
the wire pointer was measured by noting the change in the position of its shadow cast 


on a piece of sectional paper by a light at some distance behind the fixed end of the 
beam. 


Without going into the details of the measurements, it was found that the total 


torsion amounted to ‘025 radians (or 1°5 degrees) with a flexing weight of 12 lbs. and 
a deflection of about 1+ inches. 


For an aeolotropic split cylinder the accurate torsional relation is 


Wal, EF 
mE 4a “+1,) x “TOs bl 58733. 


In the present case 7 = 025/36, H =1-470,000, a= and J="11984,22489. Thus 
we have 158°1896=4#/y,+7,. The annular rings were approximately parallel to the 
plane of flexure, and if we assume y,=°"45, we have H/p,=39°44, a high value, but not 
without parallel in the case of other woods. 


For most practical purposes this method of transverse isotropy is as exact as our 
knowledge and the variability of the elastic constants of wood entitle us to use, but 
it is worth while to investigate one instance more closely with a more complete system 
of elastic constants, and to see the differences which arise according to the plane of 
aeolotropic symmetry which is taken as the plane of flexure. 


Using the figures of Wertheim and Chevandier and of Mallock the following 
systems for beechwood and pinewood were arrived at on the assumption of the ‘partial 
ellipsoidal’ elasticity of Saint-Venant. 


Taste X. Elastic Constants of Beechwood and Pinewood on 
the basis of ‘partial ellipsoidal’ elasticity*. 


Elastic Constants Beechwood Pinewood 
E 980 1113 
7 (perpendicular to annual rings) “408 “a2 
9” (parallel to annual rings) *D30 “486 
pe 172°9 45°3 
pe ta2°7 27:1 
nm: the coeflicient of ‘partial | { ; 
ellipsoidal’ elasticity ti ie 


Of course no great reliance can be laid on the accuracy of these figures, but they 
provide systems such as are likely to occur in wood and give some idea of how in 
practice the aeolotropic solution given above in § 16 would work out. 

The case dealt with for the isotropic solution at y=12° was selected for 
investigation aeolotropically—that is, the ‘projected’ problem was taken to be that 


* Moduli measured in kilogrammes per sq. millimetre. 
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concerning a non-curtate circular sector of 12° angle. The results of the calculations 


are given in the following summary : 


TaBLe XI. Aeolotropic case: Numerical results. F 
J Wa Wa 
Material Angle T/ aT ft EI for isotropic 


case with the same 
angle as for the 


o sotropic wi el 9° 0 9) aeolotropic case 
1°, Isotropic with y=}... Be hes Be yh fa eC) 846 pee 
2°. Beechwood : 

(a) with annual rings in plane of flexure aes 13° 409 6828 2837 
(6) with annual rings perpendicular to plane of flexure | 10° 313 5253 2852 


3°. Pinewood : 
(a) with annual rings in plane of flexure ee 15° 28'°6 1:9410 ‘2825 
(6) with annual rings perpendicular to plane of flexure 1 1:2226 2855 


For the purposes of comparison of these aeolotropic cases with the isotropic cases of 
a circular sector of equal angle, it is clear that the elliptical shape of the aeolotropic 
sectors is not of great importance. This is illustrated in Fig. 10. 


y = 12° 0. Projective solution, true circular sector. 
y = 15° 29". Pinewood, annual rings in plane of flexure, sector with elliptic are. 
y = 9° 18’. Pinewood, annual rings perpendicular to plane of flexure, sector with elliptic are. 


Fic. 10. Diagram to illustrate closeness of elliptic sectors to true circular sectors. 


The most noteworthy results of these calculations are to show (i) the effect of 
the direction of the annual rings in the wood on the torsion due to unbalanced shear 
couple, and (ii) the great increase in the flexural torsion when we pass from a material 
such as steel to wood. They serve to emphasise the difficulty of making applications 
of the mathematical theory to give numerical predictions regarding the behaviour of 
such a variable substance as wood. Without more precise determinations of the elastic 
constants of any individual plank we can at best give descriptive properties only and 
even then may be disappointed not to find them fully confirmed by experiment! 
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V. Summary or ANALYTICAL RESULTS OBTAINED FOR BEAMS OF CROSS-SECTIONS 
or ASYMMETRICAL FoRMS OTHER THAN THAT OF THE CIRCULAR SECTOR. 


In this section are collected the analytical discussions of the problem of flexure 
for various other asymmetrical cross-sections, which would appear to be amenable to 
arithmetical, if not to complete analytical solution. The discussion is confessedly 
incomplete in that the results obtained are chiefly in analytical forms which can convey 
no immediately apparent physical meaning. To reach the stage of physical applica- 
tion would require months of calculation and, as long as the exigencies of more 
pressing war-work do not permit of the necessary leisure, we can but look forward to 
the arithmetical completion of the work in a continuation of the present memoir. It 
seems however desirable to describe here what has been done. 


§ 20. Case of a cross-section of algebraic equation derived from assuming an 
mtegral algebraac form for the function x(x, y). An attempt to extend a method 
used by Sarnt-Venant. 

In looking for guidance in the works of Saint-Venant it is natural to expect 
that soluble cases might be found in developments of his algebraic method of 
solution™. The basis of the method is to start with the assumption that the function 
x (w, y), which arises in the expression for the longitudinal displacement wf, is an 
integral algebraic function of « and y, and to find therefrom the boundary of the cross- 
section for which this function gives a consistent solution. 


Thus, if we assume that 


x (x, y) =a, + 0,0 + by + au + 2b, xy + Cy? + 4,0" + 3b, a°y + 8¢,0y° + dy’ +... 


LBA BERS or Oe re (85), 
we find well-known relations between the coefficients which must hold if the equation 
| OX, oX 
Ox? a oy? : 


is to be satisfied. When we substitute the resulting harmonic expression for yx in 


the boundary equation (22) 
az cos (av) + yz COS (yv) =0, 
that is, in terms of the differentials dx, dy along the boundary of the cross-section 
bovlig = Ga ks OMe Re che ena (86), 
we have the differential equation of the boundary of the cross-section. 


Saint-Venant found that a large number of interesting cases, including that of 


* See the Mémoire sur la flexion (cited on p, 4), §§ 18—23. 
+ Cf. equations (14) or (20) above. 
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the elliptic cross-section, could be discussed by retaining only the harmonic a — 3ay’. 
This leads to boundaries contained in the general equation* 


m 


Cyl At 6 9F ea ee Carpe mee eset heseatees ven? (87), 


where , ©, ¢,, ¢, are, with certain restrictions, constants at our disposal. We cannot, 
however, make use of this equation for any but symmetrical cross-sections, because 
in Saint-Venant’s form of the shears yz and wz the torsion constant r was taken to be 
zero. Thereby Saint-Venant was virtually assuming symmetry of cross-section, and 


he seems to have overlooked this consideration in discussing t the geometrical forms 
m 


of contours which would be given by odd values of = in equation (87). He does 


1 
not however proceed to discuss in detail the flexure of a beam having any such 


asymmetrical cross-section. 


We have therefore to follow up Saint-Venant’s treatment with the more general 
form of equation (86) that is given in §4 above. Omitting the details of the work, 
we have the result that if we take 


) W 
ie —(B—B")a— ray + ear (n 2 (1+) m) (a@—32xy*) 00... (88), 


the boundary of the cross-section which corresponds is given by 


BoB a 2 l-m  W(1—m(1+y))(1—m) 
YY 1 —m oat & ie) nd Deraire : Ley } 
ran El! +79) (Lm) + ary Oi cent ai 


y= 0.. 


The latter equation may be written as 


where m, ¢ €{, a, b are constants which are not at our disposal until Conditions I 
and II of § 4 have been satisfied. Condition I is easily shown to be immediately 
satisfied, but Condition II presents greater difficulties, for the boundary as well as 
the integrand of the surface integrals involved depend on the constants of equation 
(90). As that equation stands, x is expressible in terms of a trinomial in y and it 
appears that in order to render the equation workable we must assume one of the 
coefficients to be small and so be able to expand in powers of this coefficient, the 
integrals of Condition II then being reducible to B-functions. This was attempted 
in two cases, (i) € assumed small, (ii) the coefficient of y* assumed small, but in each 
case the convergence of the results gave rise to doubt as to the legitimacy of the 


* This integral and also that in equation (89) are obtained from equation (86) by means of the 
integrating factor y". It may be that with other integrating factors other integral harmonics might 
afford integrable forms, but these have yet to be found. 

+ At the close of § 18 of the memoir cited above. 

{ It is the absence of the linear term it from the equation for the boundary obtained by Saint- 


Venant which forms the basis of the criticism above. 
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assumptions by failing to give appropriate values for the quantities supposed small. 
Moreover the restrictions imposed on the constants of equation (90) would be much 
more stringent than those affecting Saint-Venant’s symmetrical boundary (87) and, 
even if the analytical difficulties could be overcome, there would be correspondingly 
less interest in the solution. 


§ 21. A general method of solution by means of conjugate functions. 


The case of the circular sector cross-section with which this paper is chiefly 
concerned: might be regarded as a solution in terms of the conjugate functions given 
by the relation 

a+7B =log («+7y), 
and it is clear that the use of other sets of orthogonal curvilinear coordinates given 
by the general relation 
a+iB=f (x+y) 
might give rise to solutions which would come within the processes of analysis to an 
extent sufficient to permit of arithmetical reduction. 


With a view to simplifying the analysis we may take* 
Ww : : 
X (@, y) = — ray + Fy (1 +n) (wy? — ge’) + (a— a") y—(B,— B") &+ $b (2, y)... 


where ¢ (x, y) satisfies 


With this form of x it is seen that the boundary equation 
| -7 OG e+(1—4 )y)+B 6! di 
Y* an EL 2" ude s ly 


ox Wie 7 Ma 
Sh [re X— G2 +n)ay—(@ : ) |de=0, 


takes the simpler form — 


Ww ys ab, 
[ ary = p(t) ny) +38 | dy 3 de=... Deak (92). 


Now, by the properties of conjugate functions, 


dh _ Apia , MB _Apda _Ap2a 
oz dade OBdx dada dBoy 


0b _0pda ODOR LOR Ta. ale 
dy dady OBdy dady  eBex’ 


and along the contour given by a=constant 


da Otic, “oe 


* In order to avoid confusion of the slide constant of the earlier part of this memoir with the 
conjugate function symbol 8 we now write the former B as £,. 
8—2 
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Hence if the contour a=constant is to be a part of the boundary of the cross-section 
of the beam, we must have, from equation (92), 


da 
=” By Oat ar tata) SU (93.4) 
—s — =a 
Ox oy 
along that part of the boundary, and similarly if 8 =constant is to be a part of the 
boundary ap 
od Ox Pyle Whe aie - 
8” fa, apt areas = HY )| Shalstatee ver erates (93 b) 
ae) + Gy 


along that part of the boundary. 
The analytical problem thus becomes : 
To find a function ¢ (x, y) such that 
O° o 
~ + Fe =0 
and, therefore, equivalently by the properties of conjugate functions, 
op oo : 
ro sae 
throughout the region enclosed by a closed curve composed of parts of the contour 


: ; ; é Q 
lines given by taking a or B-to be constant, and such that ag and iy assume the 
a 


values given in equations (93@) and (936) over the parts of the boundary where 
a=constant and 8 =constant respectively. 

The algebraical difficulties that arise are 1°, to find the inverse expressions for 
x and y in terms of a and £ for substitution in equations (93a) and (936) and 2°, the 
determination of the Fourier solution for @ which will satisfy equations (93a) and 
(93b) in the way required. A failure to surmount these difficulties analytically would 
not however render the solution impossible for we can always fall back on arithmetical 
methods to solve any particularly interesting case. Indeed when it comes to the 
application of the Conditions I and IL (Equations (23) and (24) of § 4), it will be 
necessary in all but the simplest cases to use arithmetical or mechanical methods of 
quadrature to determine the value of the torsion constant 7. 


In the case of symmetry perpendicular to the plane of flexure Condition I is 
identically satisfied and Condition II becomes : 


TX [[reeuy = a L [ow — (1+) ay) dady —2(1 +) 5] 


—i]¢ (x cos (yv) —y cos (xv)) ds ...... (94). 


The integration with regard to ds must then be expressed as an integration with 
regard to a or B as the case may require, 
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§ 22. Particular cases of the solution by conjugate functions. 
(i) The Lemniscate functions*. 


These are given by 


. Cc c 
= ee a a A) hi ht ee et 8 
a+7B BE Peart reas cats (95) 
and consist of the families . 
(aot sai Pe dt engtes Ia. N, metelde Se du tilt ee wee (96) 
and est Hc Sl mi ke UE seer en hots oeerereee (oF): 


tan B 
the latter being a system of hyperbolas, each of which passes through the limiting 
points (0, +¢) or foci of the system of lemniscates given by equation (96). 


Le 


: LZ 
ar WIA 


Scale of she 


C,=centroid of lemniscate. C,=centroid of sector. #=focus of lemniscate (OF =c, OD =a=c,/2). 
Fie. 11. 


* Compare Lamé, Lecons sur les coordonnées curvilignes, Paris, 1859, p. 217 et seq., where will be 
found (p. 223) a diagram of the two orthogonal families of curves. The notation used here differs to some 
extent from that used by Lamé; it is arranged to provide a continuous system of @-values for the Fourier 
expansion. The positive half of the axis of x and origin to focus of the axis of y correspond to 6 =0, the 
axis of y beyond the focus is 8=7, the negative portion of the axis of « and the focus to the origin of the 
axis of y are both 8=2z. The lower branches of the hyperbolae through the focus are provided by values 
of 8 between 0 and z and the upper by values of 8 = 7 to 27, the two branches meeting at the focus corre- 


spond to y and 7 +y, see the diagram above. 
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When a=0 equation (96) reduces to 


(Pel 26 (ar fe parctide ss races semac eee (98). 

This is the Cartesian equation of the ordinary lemniscate 7° = 2c’ cos 26, a curve 
which takes the form of a figure of eight, the tangents at the double point (0, 0) being 
the octant lines «= + y. j 

We propose to find the solution of the flexure of a beam whose cross-section is 
one loop of the general lemniscate, the direction of flexure being in the direction of . 
the axis of « and therefore at right angles to the line of symmetry of the cross-section. 
It is easy to show that along the boundary of this single loop we must take 8 to vary 
from 0 to 27 and thus complete a circuit of the loop. 


Leb v=1—e7*cos B, u=V1 —2e~* cosB+e™, | 
then c/o =4 (u—v), (OS 2(b ED): ssavssar denne: (99), 
Ob eT ee ae aoe 
and da = a 2 (a a y’) aTy + EI (w Tt (a Wied »} ’ 


which leads to 


SB) =n tre ua (Lu) 
Kd fama, y 
We’ E aC ene A 
ST (Uy — Vy) — 2NV,) VF (Uy — Vp) (1 + Uy "*) eee ee (100), 
u, and v, being the values of wu and v when a=a,, the value of a for the required loop. 
This is the boundary condition. 
Condition I is identically satisfied. Condition II leads to . 


2pr | | ydady =" ui | | (ny? — (1 +9) w*y) daedy 


— Wy- 1 | denn (x cos (yv) — y cos (av))ds...... (101); 
see (24), p. 12, and below (387), p. 17. 
If P= (¢e’+y'-—¢)+4ea’?—ce™: ‘ 


dF /dx= 4a (a +y?+c), dF /dy = 4y (a’+y-¢), 
cos (xv) = dF /da/((dF/dx) + (dF; /dy)?)® Ate fii Bt p 
Baie? 


and similarly cos (yv) = Jape 
Thus {x cos (yv) — y cos (av)} ds = — et eds = — a eds. 
But (ds) =(de)+ (dy'= = (= : = + (“ : yt 
2 do Sa ee 2 2 20 
=f LSS agy 


ee 


ce 


oF, ds=4 on dp. 
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U 


20 2 ms Megat 
Accordingly | $u=a (cos (yr) ~y cos (av)) ds= — 4 [ Bi e~“sin BdB 
2 0 


e-* sin B 


Ce 20 
oA (Oy en Liboen if 
a), 4 ° /1 —2e-” cos B+e7 B 


Thus (101a@) becomes 


v2 | rdedy— Fy ie Uw. (np — (19) 29) drdy 2(1+n) 9 | 


Dey sin B 


eka ane. 101D). 
o V1 —2e-” cos B+ ene B ( ) 
é . re dd 
a sol : 
Now a solution of ar dp? =) 
is p= (A, sinh na+ B, cosh na) sin n(B—7).........40. (102q). 


Inside the loop a=0, a ranges from 0 to + at the focus (0, ¢), but w and 
therefore ¢ must be finite; it follows therefore that A,+B, must be zero and our 
solution will take the form 


PaE (dye “sin m (Bom) bce hectic ease (102), 
(3) EPS (nA aan (Bor) ee a, (103). 
BaF oe kaiiciye fiom (ho0v and (103) ty 

 (nd,e-sin n(B—n)) =47e°/(B) + e407 (B), 
where f(6)=Vur—af(14+%7), ¥(8)= (4 —4—29y) VE maw) (1414-9). 
We 


If viz = ave +444 a EI Bisa 


£8) sin n (B — 1) dB, 


nN 


Coe 
nut + sin 27 


A wee ie c W(B) sin n (B—7) dB. 


Nar ar -+sin na sin 27 


it follows that A,/ = 


These integrals must be obtained by quadratures and the values of 4,’ and 4,” then 
inserted in the f-integral term of the torsion equation (101), where again for the 
general case quadrature is needful. The proper values for can only be determined 
when we have found the expressions /(8) and (8) for the particular case. 


We will illustrate the method on the relatively simple problem of the torsion 
flexure of a beam, the section of which is the loop of an ordinary lemniscate. In 
this case a,=0 and 

e-%sin B 


ie Pa= OM eeat aoe 1B | “danocos $BdB, 


2e~ % cos B + e - 2% 


Now in this case 
%=2sin?4B, u,=2sin $8, 
and we find 
a=csini(r—8)J2sin3B, y=c cost (r— 8) /2 sin $B....... (104) 
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over the contour a=0, and accordingly 


—47c (sin B+ cos $8) 


ae 
a a 
aie 
2 
Il 
oO 

ll 


Ty 2 sin FB {cos }(38 — 7) +($+)(sin 4 (8-7) —sin}(5B—7))}...... (104b). 
Thus f (8) =sin B+cos 4B, 


(8) =V2sin $8 {cos } (38 —m) + (4+) (sin 4 (B—7) —sin 4 (58 —7))} 


A suitable form of ¢ is 


p=4tre (e-* sin B+ 2e-4+ cos $8) 


w We’ ; 
+> SA," +(4 t+) A,'’) 1BI*? sin 2 (B— n)| Bae (105). 


This will satisfy (1046) provided we take 


Ayre F /2 sin 4B cos 4 (38 —7) sinn (8—7) dB 
; ..(106). 
Ae a '/2gin sin $8(sin 4 (@—7) —sin4(58 —7)) sinn (B—7) a) 
These may be shown to be 
Hie = /2 cos 18’ sin 28’ sin nf’ dp’ 
senate (107) 


i Moe 
Fi (2 cos 4’): sin 38’ sin nB’dp’ 


nt 


These integrals were found by quadratures, using one of Sheppard’s formulae (Bio- 
metrika, Vol. I, p. 276, Case (i) (c)), the subject of integration being tabled to every 5°. 
After considerable labour there resulted : 


; ae Aa aC kp=t4g + EA | y= SA 
ee 
1 | —1-125,0142 | +1:500,0012 | —-333,3333 | —-093,7535 | +-375,0003 
2 | + -031,2632 | + -187,4978 | +-133,3333 | 4+-031,2530 | + -046,8745 
3 | — -013,0284 | — -020,8406 | —-085,7143 | —-005,8622 | —-005,2102 
4 | + -006,8490 | + :005,8517 | +-063,4921 | +-002,4437 | + -001,4629 
5 | — -004,0996 | — -002,3637 | —-050,5051 | —-001,3204 | — -000,5909 
6 | + -002,6992 | + -001,1255 | +-041,9580 | +-000,8155 | +-000,2814 
7 | — -001,8626 | — -000,6602 | —-035,8974 | —-000,5482 | —-000,1650 


Now ¢a=0 =4rc’ (sin B+2 cos $8) +3 (An bd 60) ea) ae UG sin 2 (B — 7) 


4EI 


and le fa=0 cos $BdB =4rc° ($ SETS et 


; iy, 2 
TET (40! +(£+7) A,’ yt 
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Equation (1016) for the torsion now becomes 
T 2 [[ydeay —4e' (B+ 2n)| 
W 
se gil | {tw -( +) wy; daedy —2 (1+) Ly 


os f my a ie n(—1)" 
$03 4(4," + (b+9) 49") \ ba (108). 


But for the common lemniscate 


[[y diady = $c (L4 +), [|r dy =t37¢ lie y dady =", 1c’, 


thus we have 


Wef{ w/74 @ tr 29\ 
( 33) 


el ee i ny a br 3 yal : 2 652 
7(—47) EY \alag~ 4/72 (a +°381,1119 —(4+4 7) -472,65 at 
We 
t (~°7853,9816) =5,, {02,4468 + 121,31947}, 
We 


7 = — 7 {08,1154 + *154,46879} 2 .ceecseseeeenees (109) 


the breadth of the lemniscate. If 7 =4 


Ox 02951. 


T= 


ee, 
Thus if the lengthy analysis be without slip, there is a small negative torsion, the first 
we have come across. 

We can now complete the solution. We will write 


We ( ) 
T= — ET €c3 7p), 
and we have 


(x, y)= - moe: «sinB +e-4* cos $8) (e, +e) 


3 {(fA," + 4A," +49," e-" sin n (8 — mY] 
1 


=~ Fp7 [(dent sin B eH 00848) 6 +en) — 5 { (hu taaln)e-msin n(B—m)}] 


the numerical values of x,, x,’ being given on p. 64. Hence by equations (20) 


and (91) 


ay? —% ges 
ces aN (es) . J 
ies 2) 4 3+ myn) e-*sinn ate — Bat ay +7: 
: 1 A 


u and v will be as in equation (20), p. 12. 
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We have now to determine the fixing of the clamped end; let us suppose the 
element at the centroid fixed. Then, by p. 42, a,=,’=0 and 


Wree 
y=—ayy, 7 =(aten), 
Wiel ies 
a2 py ens — vy +e = 0. 
“y oe VY a cty LIA Vos 
Hence u(0, y= — Tye + pr len ye t gle — 62] + Bex 
: . W 
= py (le 42) +2(B,—1I4+ 407 oy) 
= iat Ls 
on SB ytd py 
gives the linear term in the deflection as soon as f, is found. 


We have now to find o along the axis of symmetry. Clearly 5 =0, when x=0, 


OG ees = , when «=0. 


da Ce 


But e-*=(—y)/e for y<c, =(y'—©)/e for y>e, and at the centroid the former 
must be used for y=47ce. Further 8=0 for «=0, from the origin to the focus, and 
=7 beyond the focus, 1.e. y=c. Thus between origin and focus 


dw We [eae pees ¥ _ a—a? 
Te = ET | eten) (2— Y\+4n% oe 


3 3 2n(— Ihe O08 (4/0) | z 
1 


Further, by equation (20) along «= 0, 


du Ww 
bn Tt plang’ + 4—322)+ B.. 


Hence adding 


x2 . dw du iy We y y’ 
(F) <2 ~ det ps PORE Ne (a+ én) = pe cm 


- £3 2n(— 1)te= O08 (4, + myn) | per (110), 
which enables us to get the shearing stress from, the origin to the focus. Beyond 
the focus B=7, thus cos 8 changes sign, and c’—7’= —c’e-4, thus the term }e-¢sin B 
in w remains the same as before in dw/dx, but the last term changes and we have 


we We . 
(=).-0 = aT | + €,7) Z 49% Y uP Ys > Qne~("-1)4 (ie, + ky ) | <t CL pO) 
lyse Bi 


ON THE TORSION PRODUCED BY FLEXURE 67 


At the focus a= and the summation reduces to a single term the same for both ; 
thus we have 


xz’ We 
(F)-0 = ay leat 2m t+ (e+ 1+ 2x) 7] 


pee told 
=" ay | — °184,8916 + 1:904,4693m]........-.00., (111) 
We 1 Wa’ a 
= py X 29178, ifn =4, a hati 14586, if a=J/2e. 


At the centroid we have 
N72 2 = 9—(n—1)a ( _.]\n 
B,= We E NS ie dds ety ese {0 


Eicon 16 einen 2 
and e-*= (8? — 9 )/e = 1 — yaa =°3831,4973. 
x2 We a ane~(m—1)a(—~])n 
But =) oye | yi iveadt 
u (= ile wie (ten) 4ety 7, — = (kn + Kn’) 2 
Wy  _ 
oc =43y Te T+ Bex 5 
or, the curve of deflection is Jel 
W xz \ 
u(O, 4 —- (4/2 —14") 42 (= oe Onan sete Tain eriee cas 1 he 
( y) = Rp 6 ) in lh=5 ( ) 
as it should be. It only remains to find B, and s,. 
Numerically we have 
n Zane—(—1) a ( — 1)" Ky X krne—@-1)@(—1)” | Ky’ x$rne Na ( —1)" 
1 — 1°570,7963 + °147,2677 —-589,0491 
2 + 1°:203,7003 + '037,6192 + 056,4228 
3 *691,7963 + °004,0554 + °003,6044 
4 + +3! 53, 4154 + -000,8636 + °000,5170 
5 - 169,2637 + 000, 2235 + :000,1000 
6 + °077,8238 + °000,0635 + °000,0219 
7 | . — -036,3589 + -000,0199 + -000,0060 
These values lead us to 
; We’ . 
B,= spy [308,425 19 —°190,1128 + °528,37707] 
= C 
arln LOD IMAG CSHB ST yl ee tetetlsants cl bstvaess (113) 
ie a 
—> x ‘009,544, 
Sie) Sal 
if 7 =4 and a= 2c be the breadth of the lemniscate loop. 
We 
Further, tay (1 266,546777)—e LB PGE GO | ace cen wien ses (114) 
We? 
z= 449, 
EIT x 06 


if =4 and a=V2c be the breadth of the lemniscate. 
9—2 
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Thus the deflection will be 


ee Me 3 Gj) -3 (3) +3 sa 06449 (| Eas ae © (115). 


The linear term is only about 1th of the linear term in a circular sector of angle 90°. 


This is due to the change in the torsion from a sensible positive to a slight negative 
value. 

It remains to investigate the vertical shear along the central radius. In working 
this out attention must be paid to the change in form in the shear after passing the 


focus. The following table gives the chief results : 


Shears parallel to load along Aais of Symmetry. 


n=% 
—~ {Wat 
yle y/a za | EI = [ee =I cae 
tz | ar xz/(mean shear) — 
0 0 0 0 —1 
al ‘OT 07 019,541 + -038,289y — 0100 — 1135 
2 1414 037,932 + 087,573 — 0160 — 1-216 
| “4 2828 067,834 + °222,562 — °0122 — 1-164 
6 “4243 086,317 + °410,126y + 0162 SP 
43 DHD4 093,833 + °633,273n + 0645 — ‘131 
1 ‘7071 092,196 + *952,235y + 1459 ,+ 966 
115 8132 — 084,637 + 1:214,4149 + °2190 + 1/951 
1°30 9192 — ‘071,986 + 1:509,133% + 3053 +3114 
J2 | 1-0 — 060,879 + 1°754,884y | +3778 + 4-091 


The last column is plotted on the diagram, p. 61. It is clear that from the origin 
onwards we have a small negative shear and ultimately a moderately large positive 
shear, reaching a maximum at the butt end of the section. The form of the dis- 
tribution differs widely from that of the 90° circular sector. We have not succeeded 
in discovering any physical reason for the torsion being negative in the case of the 
lemniscate and positive in the case of the circular sector. On the other hand, we 
have not discovered after much investigation any slip in the analysis. We publish the 
investigation in the hope that a fresher eye may discover our slip or verify the analysis, 
or again find other sections with negative torsion. Should such sections actually 
exist, the quick transition in sign with changes of not very marked character in the 
asymmetry might suggest that the edge of attack in a propeller-blade may be of a 
rather unstable nature. Possibly when time is available for the sections suggested 
below to be worked out numerically, light will be thrown on this interesting 


question. 
(il) The coaxial circle functions*. 


lhe equation aeelee x +t (y +e) 
-e2+2(y—c) 


1 \ . . nuts 5 i Fy - ro 7° 
* “Systeme cylindrique bicirculaire’ of Lamé: see his Lecons sur les coordonnées eurvilignes, p. 199, 
Paris, 1859. 
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gives rise to the two families of circles 

e+ (ye coth a)’ =¢ cosech? @ ......6..ccccceseeeee CLLZ); 

(x —c cot B)?+y?=C' cosec? Bo o....cscccevssseceees (118), 
the a-family being a system of coaxial circles having real limiting points at (0, +c) 
and the axis of « as radical axis, and the B-family being the orthogonal system of 
coaxial circles. 

The interesting case from the present point of view is the boundary given by 
taking two circles of the a-family, one being wholly inside the other*. The beam 
then is a hollow shaft with an eccentric cavity, and thus possesses the requisite 
uni-axial symmetry. 


The expressions we require for substitution are : 


est fa) 
| i or tig Ba eaten (119) 
csinh a 
4 cosh a —cos | Be 
| eee 
0x Cc 
da _ 1—coshacosB | (120), 


oy ie mAh ibis coneah tebe wea 
(=) (=) = (cosh a —cos f)? | 
Ba NB : 
se g C 


and it is easily seen that a complete circuit of any circle of the a-family is given by 
taking 6 to vary from —7 to +7. 
Equation (93) becomes for this case 
op orc? sinh* a sin 8 We’ (1+) sinh asin’ 8 —y sinh’ a sin B 


‘ (cosha—cos8) EI (cosh a — cos £)* 
Dith ks eam (121), 
and the form of the solution will bet ; 
ss 2 A,, cosh 7 (a,—a) +B, cosh n (a, —a) Sos (122), 


a sinh 7 (a, — a,) 
where the constants A, and B, are determined by Fourier’s method from the equations. 
obtained by substituting a, and a, for a in equation (121). 


Condition I is again identically satisfied. Condition II becomes 
+ [yrdedy = ue ae +) wy) dady —2(1+7) | 
tye] " SNS ECP ake T 55). 


cosh a—cos 8 
* This problem is suggested as a soluble one by Love (Z/asticity, Second Edition, § 232 (e), p. 325). 
The corresponding torsion problem has been investigated by H. M. Macdonald, Proc. Camb. Phil. Soc., 
Vol. vill, p. 62, 1893. ; 
+t Lamé, loc. cit., p. 184. 
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In the B-integral a in ¢ and cosh a must be given the two selected values and the 
integrations round the two circles taken in opposite senses. As in the case of the 
-lemniscate ¢ will contain series terms with both 7 and W/EZ as coefficients. 


To get the solution for a beam resembling in cross-section the blade of a propeller 
we might use the curvilinear quadrilateral enclosed between the circular ares given by 
a,, a, B,, —f,, but in the sector of an annulus we have such a similar figure that 
there would be little but the pleasure of confirmation in following out this case. 

(iii) The confocal cone functions™. 


The appropriate systems for our problem are given by 


a+78 =sinh"' ° (ce + iy) Sek ece nts Lee. (124), 


c 
and consist of the family of ellipses with foci at (0, +c) 
oo y 


a CORD. == (Sud ieis.0 a efeleteinie ins ep oleinie tel aieteisiarane (1 25), 
and the family of hyperbolas with the same foci 
adh os Sabet 
cos B aro Bo ee (126). 
The Cartesian coordinates of the system are 
«=c sinh a cos B) (127) 


y =c cosh asin B) 
and it is seen that the part of the axis of y lying between the foci may be taken as 
the ellipse a=0, and the part lying outside of the foci as the hyperbola 8=4z, and 
the positive axis of « as the hyperbola 8=0 and the negative axis of x as the hyper- 
bola B=7z. 

From the point of view of the flexure problem, and especially in relation to the 
problem of the propeller-blade, it would appear that the most interesting case would 
be the boundary enclosed between an ellipse a, and one limb of a hyperbola given 
by taking B=5-B, where 8, has an appropriately small value so that the eccen- 


tricity of the hyperbola will be fairly large. As in all cases, however, where composite 
boundaries are made up from systems of orthotomic curves, there is the disadvantage 
of the right-angled corners which do not appear in the section of a propeller-blade. 

For this case, where —a,<a<a, and 3 —B,<B< set B,, the boundary conditions 
derived from equation (93) are easily found to be 


0 ‘ ; ; : 
as = 27’ cosh’ a sin B cos 8+ ie ((1 +7) sinh’ a cosh a cos’ B—y cosh’ a sin’ B cos f) 


* ‘Cylindres homofocaux’ of Lamé, loc. cit., p. 195. 
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and 
od ari fo} We anlietie A } hie 
36 al 27’ sinh a cosh a sin® B — Er ((1 +7) sinh’ a sin 8 cos’ B — y sinh a cosh’ a sin’ 8) 
Serr ie Pera ee ae (129), 
and the solution for ¢ will be of the form 
ee n arent) Fe n 
=A, sinh a (B —- _ 6.) cos — (a—a,) +B, sinh 3. (a —a,) cos a (8 _ : - p,) 
DA tes aes (130), 


where the constants A, are derived by Fourier’s method from equation (129) for 


B =5+ 8, and the constants B,, from equation (128) for a=a,. 
Turning to (94) the [e (x cos (yv) — y cos (av)) ds = —Leé [sin 2B¢.,¢8 along an 
_a=a, boundary, and to +4c’ |sinh 2ad¢g,da along a 8, boundary. Thus the equation 


for the torsion is provided as soon as A, and B, have been determined. 

It is evident that this solution can be completely set forth ready for arithmetical 
calculation without methods of mechanical integration. 

(iv) The rosette functions. 


We shall lastly consider the system of conjugate families given by 


When polar coordinates (x =r sin 6, y=7 cos 6)* are used the equations of the 
two families are found to be 


a =— cos mé | 
Se ht La ier a ee (132), 
B=— sin m6 | 
and we easily deduce that 
roe as io et 
w=a(a'+ fp’) sin & tan”? | 
aR Cire eA (133), 


: 
1 

y =a (a? +B’) ™ cos & tan“ e | 
and that 


CL qinrt 


a te da ee ~~ ma” 
ox oy 


= -= (a? + By (eB cos (F tan™ 4 +asin e tan” e).. .(184). 


a 


sin (m+ 1) 0 


* An unusual form of polar coordinates has been adopted in order to preserve uniformity with the 
notation given in the earlier portion of the paper. 
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The typical curve. of the a-family is a rosette composed ‘of petals radiating from 


1 
the origin and enclosed within a circle of radius aa ™, which is touched by each petal. 
Let us take the rosette a=i and let us consider the single petal which results from 


the principal values of tan~ 8 being used, that is for ~5 <tan"*B< + This petal 


2° 
is a closed curve symmetrical about the axis of y and having a node of angle = at the 


origin. With a suitable choice of m we thus have a curve which would resemble fairly 

closely the cross-section of a propeller-blade, and it is with this in view that we proceed 

to give a formal solution of the flexure problem of a beam having such a cross-section. 
The boundary equation (93) becomes after reduction 


Gf) -- Rae +B)" * cos ( tan”! B)+ ula (1 EY oe 


x (a +7) sin’ & tan™ B) — 7 cos’ & tan” a}) | ( cos = * tan” 8) +sin a ) tan ‘B)) 


and this is to be satisfied for the whole of the range —-« <B<o, a range which 
renders a Fourier serves solution inadequate. 


Consider the integral 


ea ua 1) 
a=2{_ |, “P| (w(- — — — 7, e081 (v=) dud ...(136). 


The integrand is ames “ if F he: is finite for —« <v<o, and it will 
be assumed for the present that the integral is convergent* for the range 0<a<1, 
—«o<B<o. Also it satisfies the differential equation 


od od 
dai * age 
and we have 
(=) 20. ” Fv) 0s u (v—B) dudv (137) 
Bai) | EO CORR BND ae) tHe 


=F'(B), by Fourier’s integral theorem. 
If, therefore, for F' (v) in equation (136) we substitute the expression for 


aa) 


0a 


—u(a—1) 
* It does not appear to have been noticed by Fourier that the integral | z ea cos u(v — 8) du is 


divergent for any range of integration which includes the point w=0. Compare the Théorie analytique 

de la Chaleur, § 421, Darboux’s edition of the Wuvres de Fourier, Vol. 1., p. 507, Paris, 1888, where the 

integral 

i °) 2) dp “5 
daf (a) | Mi cos (pa — pa) (e?Y + e-PY) 

is considered without any mention of convergence; this has probably been often remarked before. The 
i) oo —u 

convergency of the integral (123) above is obtained by the insertion of the constant | | F (v) — dudv 
-x2 J) 


which is also infinite. 
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from equation (135) we obtain 


¢ (a, B) =i [. (“ (a )) (S ~ et cos (v -B)) dudv ...(138) 


T 


as the formal solution of our flexure problem, when combined with equation (91). 
The expression for ¢ (a, 8) takes a simpler form when we put v=tan w, 8 =tan ¥. 
The value of 


| {x cos (yv) —y cos (av)} dds 


in the torsion equation (94) reduces to the simple expression 


Ti ay 2 (+2 1 
m al Agen) Bg _ 5 [4.4 aaa: 
il +B’) m 2 (1 + By” 


a? +m fo 0) +0 ‘ mel 
St | | | F(v) sin u(v—B)(1+ 6’) “dBdudv, 
0 -2 


eS 2) 


which completes the formal solution, but is an expression hardly likely to prove of value. 

The portion of the integral (136) which concerns w depends essentially on the 
function log ((a—1)’+(v—£)’), and it is suggested that a variant formal solution 
would be obtained by taking 


$ (a, 8) =1{ Fe) foe (ea ty cee Vier eee (139). 


Assuming convergence, we can easily show that this integral satisfies the 
differential equation 


and that 


so that as before we can take as solution 
Ted le (e{e, »)) 
b(a p= —{ (Eo 


The question of the convergence of the integrals used in this discussion clearly 


log {(a—1)' + (v— B)"} dv seeped (141). 


a=1 


requires further investigation. At first sight it might appear that the integral in 
(139) is divergent, seeing that logarithmic term contains (v—£)’, but it is to be 
remembered that /’(v) is in our case an odd function of v, so that it is not decided 
that 

| F (v) log {((a-—1) + (v—B)’} dv = Lt | F (v) log {(a — 1) + (v— B)’} dv 


is infinite. There is further to be decided the uniform convergence of the integrals 
when B=+0. These points are left over for future discussion. 

The crude formal results given in the latter part of this section may prove on 
further investigation to be unworkable, or at least unprofitable, and are published 
here only by way of indicating the lines along which it is intended to continue the 


work when opportunity offers. 
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CoNCLUSION. 


To the engineer and especially to the designer of aeroplanes the chief interest 
of this work will lie in its application to the problem of the change of the ‘angle of 
attack’ that arises in the blades of a propeller rotating at high speed, and it may be 
well to state here briefly the results that seem to affect that problem. It is first of 
all to be emphasized that the solutions given in the memoir are solutions of problems 
of much less generality than those of an actual propeller-blade and that the results 
can therefore be regarded only as first approximations to what is wanted. Thus the 
beams considered are uniform in cross-section throughout their length and that cross- 
section is in every case symmetrical about an axis which is perpendicular to the plane 
of flexure. In beams of cross-sections having no axis of symmetry the torsion effect 
would presumably be of greater importance*. Moreover the loading is not continuous 
as it is in an actual propeller-but consists of a single load acting at the end of the 
beam. With these reservations we can state the result : 


For a beam of isotropic material and of uniform uni-symmetrical cross-section of 
shape resembling the cross-section of a propeller-blade, flexed by a load of weight 
W applied at one end in a direction perpendicular to the axis of symmetry, the 
torsion per unit length due to the unbalanced shear-couple may be taken as approxi- 


mately + equal to ‘28 gue when # is the Young’s modulus of the material, a is the 


ET’ 
breadth of the beam and J is the second moment of the cross-section about the axis 
of symmetry. The direction of the torsion is such that the thicker side of the beam 
twists upwards with a downwards flexing force. 


In beams of cross-section resembling the cross-section of a propeller-blade the 
torsion due to anticlastic curvature is in the opposite sense to the shear torsion and 
roughly about half its magnitude. 


If the material of the beam be aeolotropic we have seen that the torsion due to. 
the unbalanced shear-couple may be increased from twice to six times the isotropic 
value. This increase is largely determined by the ratio H/p, where F is the longi- 
tudinal stretch-modulus and p is the transverse slide-modulus parallel to the direction 
of loading, but it is far from being directly proportional to this ratio. 


When, however, wood is the material, each particular beam would appear to have 
its own amount of torsion, the asymmetry of the elastic properties of the material 
being apparently more effective than the asymmetry of shape. To the designer of 
aeroplane propellers this should be of great interest for, as mentioned abovet after 


* The amount of torsion will, of course, be a maximum when the load is perpendicular to the axis of 
symmetry. If it make an angle wy with it, only the component W sin w produces torsion, 

+ For the calculations on which this estimate is based see $§ 10 and 13. 

+ irs 19; 


ON THE TORSION PRODUCED BY FLEXURE 75 


the description of an experiment which brought out the importance of this asymmetry 
of texture, propellers might be constructed of layers of wood selected and arranged 
in such a way as to minimise the torsion which arises from flexure. 


Apart from the propeller-blade approximation the most interesting case theoreti- 
cally is that of the split pipe in which the torsion effect is so pronounced. We do not 
know of any instance in practical engineering in which this is of much importance 
although such may well exist, but it is conceivable that some practical value could 
be found for the phenomenon in such a pipe, for example as an economical form of 
cantilever weighing machine for heavy loads or in circumstances in which the ordinary 
methods would be inconvenient. It might also be used for light loads as a letter- 
weighing physical toy. Such value arises from the fact that.the angle of torsion is 


‘always proportional to the suspended load and is very large if the wall of the pipe 
be thin. 
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